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I. INTRODUCTORY CONCEPTS AND TOOLS 


la. Scattering Matrices and' Scattering Cross-Sections 

In scattering experiments of any kind, the fundamental 
observable quantity is the scattering cross section; and the 
theoretical object of prime importance associated with it is the 
scattering matrix (or S-matrix) . We begin by giving formal 
definitions of the S-matrix applicable to general scattering 
processes and in later sections ' adapt them to molecular colli- 
sion problems. The S-matrix will be shown to be directly re- 
lated to the scattering cross section through an auxiliary 
object, the transition matrix (or T-matrix) ; and a key objective 
of the theoretical program of collision processes is the calcu- 
lation of S- or T-matrices. 

S- and T-matrices are quantum mechanical objects and 
they will be introduced as such in the present subsection. How- 
ever, under certain conditions, such as in the treatment of 
heavy-particle dynamics in molecular collisions they can be 
adapted to semiclassical or even classical evaluations . In the 
present lecture series we shall fochs on the semiclassical 
evaluation of the' S-matrix (called the "semiclassical S-matrix"), 
with special applications to the case of electronically non- 
adiabatic transitions. Before discussing this explicitly, we 
shall review some basic ideas associated with the S-matrix in 
general . 

,A collision system^c|n be thought of as being in the 
states (incoming) and 'P (outgoing) before and after the 

collision process respectively? and the collision can be visu- 
alized as a "blackbox" mathematically described by a transformer 
operator S (the scattering operator) such that 

. (la.l) 

Suppose 4'^'^ is specified by a given free state | i> (this state 
is often-^chosen to be an eigenstate of the unperturbed Hamilton- 
ian) . can be expanded in terms of the complete set of 

which I i> is a member 

= I c. ! j> . (la. 2) 

j ^ 

The transition probability amplitude from | i> to a particular 
j f > is then .given by 

= <f]¥°'^'^> = <f|s|'l'^^> = <fls|i> = . 

(la. 3) 



Hence, it is seen that Sf^, which is a matrix element of the 

scattering operator S in the [ i> representation, gives directly 

the transition probability amplitude from the state | i> to the 

state I f > . We note that the conservation of probability for 
,„xn . ' ...out . , , . 

'i' and '¥ requires that 


oi-outuout^ 



(la. 4) 


implying that 

i" i 
S S = SS =1 , 


(la. 5) 


which is a statement of the unitarity of the S-matrix. 

The definition of S can also be formulated in the time- 
dependent formalism, if we set 

= 'F{t=-«) 

,pOut ^ . (la. 6) 

Then S can be defined by 

¥(+“) = S'?{-“) , (la. 7) 

and is also expressible as 

S = lim U(t,t ) , (la. 8) 

t -> .5- oo 

where U is the time evolution-operator for the interaction pic- 
ture wave function. 

In order to calculate the scattering cross section from 
the S-matrix, we first evaluate the transition probability per 
unit time since the cross section is defined as 


CJ 


i^-f 



(la. 9) 


where F . is the incident flux of particles for the state [ i> . 
Now S, as defined earlier, is the transition prob^ility ampli- 
tude between infinite past and infinite future times . To find 
the probability per unit time, we must express it as the limit 
of a time-dependent quantity S£j_(t) for i?^f (transition between 
different states) such that 


S-. = lim S . (t) . 
fi fi 

t-> +“ 


(la. 10) 



The transition probability per unit time can then be de- 

fined as 


W. ^ 


lim 

■t>0O 


|Sfi(t) 


(la. 11) 


Since we are interested in transitions between different states 
I i> and [ f > , it is reasonable to extract from S a term which is 
zero when i and f are indeed different/ namely AlsO/ since 

energy conservation holds between infinite past and future 
times, we can further separate out an energy conservation 
factor 6 (E^-E^) in the remainder term. Hence we write 


S^. = 5^. - 2iri6 (E.-E-)T,. . (la. 12) 
fl fx X f fx 


The factor 2iri is introduced for convenience; and T, whose 
matrix elements are T^^ = <f|Tji>, is called the transition or 
T-matrix. Written as in (la. 12) with an explicit energy depen 
dence, Sf^ (for ij^f) can now be expressed as 


fi 


- 27 ri 6 (E.-E^)T^. 

X f ' fx 


S .(t) 
t^-“ fi' 


(la. 13) 


This identifies S^. (t) as 
fx 

S_. (t) = - T 
fx ^ fx 


t/2 


dt' 


-t/2 


i(E^-E^)t'/h 

e 


sin{ (E , -E ) t/2h} 

_ 4 - m X f 

-h ^ ^fi (E^-E^)t/2Ti 


(la. 14) 


which represents a transition probability amplitude for i-^f 
over a finite period t of observation. The transition prob- 
ability per unit time, from (la.ll) , is then given by 


^i->f 


It,,!^.. sin^{ (E.-E^)t/2H} 
fx Ixm X f 






Ei-Ef.2 


(la. 15) 


2fi 


which is equivalent to 


6(E.-E^) |t^. I 
x->-f -h X f ' f X ' 


f^i . (la. 16) 


In general, in any experimental situation, there will be a group 
of final states | f > satisfying . energy conservation, and the 



scattering cross-section observed will correspond to the sum 
over these final states. From (la. 9) , 




^ I W 

Fi ^ ■ 


i^f 


(la. 17) 


The sum will go over to JdEjP (E^) if we are dealing with con- 
tinuum states in the product channels; p (E^) is the density of 
states - 

Before leaving our formal discussion of the S-matrix 
we will point out a formal relation between it and the Green's 
function to set the stage for the introduction of the semi- 
classical S-matrix in Section Id. 

The operator (E-H+in)”^ (n being infinitesimal) is 
called the retarded Green's function operator and is directly 
related to the T-matrix, and hence the S-matrix: 

g'^(E) = Gg(E) +Gq(E)TGq(E) , (la. 18) 

where Gq (E) = (E-Hq+Iti)”^ is the imperturbed retarded Green's 
function operator. The Fourier transform G'^(t) of G"*'(E) is 
given by 

dE e^^'^ (E-H+iii) • 

(la. 19) 


G'^(t) = 


1 lim 


2f!fi n-K) 


By actually evaluating the integral in (Ia.l9) , G^(E) can be 
written as the inverse Fourier transform of G"*" (t) : 


g'^(E) 


(±/n) 


r 

dt exp (iEt/R) exp (-iHt/K) . 

•'o 


(la. 20) 


(la. 18) and (la. 20) will be the starting point for our formula- 
tion of the semiclassical S-matrix in Section Id. 


lb. Partitioning of the Molecular Collision Hamiltonian 

Beginning with this subsection we will restrict our 
attention to molecular collision problems, and the general 
formulation of the S-matrix given earlier will be adapted 
accordingly. To prepare ourselves for these special formula- 
tions of the S-matrix we will start with a discussion of the 
molecular collision Hamiltonian. 

Separation of Electronic and Nuclear Motion and the 
Born-Oppenheimer Approximation . Consider a molecular col- 
lision system with n electrons and N nuclei, with electronic co- 
ordinates denoted by rj_ and nuclear coordinates by Rj_. The 
total Hamiltonian can be written as 


where Hgi (R/?) stands for that part of H besides the nuclear 
kinetic energy operator. The first-step to simplify the Hamil- 
tonian is to separate out the center of mass motion of the whole 
system. The dynamics of the- system is then described by an 
arbitrary set of (3N+3n-3) independent internal coordinates. 
Several choices of these can be made to suit particular problems, 
but here we will not be concerned with the specific choices . In 
general one attempts to make coordinate transformations which 
would leave the total kinetic energy operator in the form 


T 


„2 ^ ^2 R. 

r 1 

2 M S 2 4 y. 

1 1 


I + (m.p.) (Ib.2) 
^e i i 


where (m.p.) represents mass polarization terms (proportional to 
Vr -Vr and ) which can be ignored. In (Ib.2), M stands 

i j i j ^ 

for the 'total mass' corresponding to some C.M. coordinate S; 
yj_ and are the reduced nuclear and electronic masses corres- 
ponding to the nuclear internal coordinates R^^ and electronic 
internal coordinates r^ respectively. The total Hamiltonian . can 
be written in terms of the transformed coordinates as 


H 




2M 


, N=1 

^ I 

2 i=l 




+ V(R,r) . 

i=l- ^i 


(Ib.3) 

V(R,^ may be taken to include spin-orbit interactions. The 
Hamiltonian in terms of the 3N+3n-3 internal coordinates be- 
comes 


H = 





+ V(R,r) 




H (R,t) 
el 


(Ib.4) 


where 



(Ib.4) looks identical to (Ib.l) , except now ^ and 'r are under- 
stood to be the internal, coordinates with the C.M. motion 
separated out. Hamiltonians will be written in this sense for 
all our later discussions . 

Separation of electronic and nuclear motions ultimately 



depends on the great disparity between electronic and nuclear 
masses. Since the forces acting on them are of comparable 
orders of magnitude, nuclear motions tend to be much more slug- 
gish than electronic motions and electrons can be assumed to ad- 
just to the nuclear motions adiabatically . This consideration 
is the basis of the Born-Oppenheimer approximation, which in 
general means the asstimption of the existence of some basis set 
of electronic wave functions (r";^ which depend parametrically 
on the nuclear positions and for any fixed ^ is orthonormal 
and complete: 


d^r^*(r;^^^(r;R) = fixed R , 

(Ib.5) 

i “ Z fixed R , (Ib.6) 


such that the total wave function ^ (R,r) of the collision sys- 
tem can be expanded as 


i|j{r,R) = 1 x^{R)4>^(x;R) . (Ib.7) 

I 

Adiabatic vs Diabatic Representation - Nonadiabatic 
Couplings and Transitions . The nuclear wave function 
(R) describes the motion of the nuclei on the potential s\ir- 
face associated with electronic state £. We will see below what 
we mean by these surfaces when different basis sets are 
chosen. Specifically, we will see hov? the different choices will 
lead to the adiabatic and diabatic representations. The func- 
tion Xji(R) will asymptotically give the amplitude for the 

probability that the system ends up in the electronic state 
and thus directly gives the S-matrix for this transition (whose 
initial state is specified by the problem). Hence, our objec- 
tive of obtaining the S-matrix reduces, quantum meclianically , 
to the solution of the Schrodinger equation for (R) • 

The time-independent Schrodinger equation to be solved is 


Hijj(r’,R) = E\p(r,R) 


(Ib.8) 


where H is given by (Ib.4) . S;ibstitution of (Ib.7) into (Ib.8) 
results in the coupled set of equations for x^ • 


(T„ + T'.'. + U. . - E)x- = 
R 11 11 '^i 


- z 




(Tl . + TV . + U. .)X- 

1] i: ID ^0 


(Ib.9) 



where 


u. . 

“ <9. 


N=1 

t! . 
13 

= 1 
k=l 


N-1 

^11 

ij 

= 1 
k=l 


<®. |h , U .> 

^i' el 

2 >, 


-■n 




2y, 


<^ . V (fi .> • V 


<!>.> . 

X Rj^ : 


(lb. 10) 
(Ib.ll) 

(Ib.l2) 


Brackets denote integration over electronic coordinates r only. 
The diagonal elements Uj_j are the - effective potential energy 
surfaces mentioned earlier that govern nuclear motion. The 
diagonal terms T^j^ , which are nonadiabatic corrections to the 
potential energy surfaces, are usually small and can be ignored. 
From (lb. 9). it is seen that motion on the different surfaces 
are coupled by the off-diagonal elements Tj^ . , T^j and Uj_ ^ . T| j 
and T^ j are usually called nonadiabatic couplings , of whxch 
T|j are the dominant ones (dependent only on the first derivative 
of (j) with respect to nuclear coordinates) and are also nuclear 
velocity dependent (since it is proportional to operating on 
nuclear wave functions) . 

The adiabatic representation is defined by the requirement 
that the potential matrix U be diagonal for a!ll R- This repre- 
sentation employs for the electronic basis functions the eigen- 
functions of H , : 
el 


the superscript a on the wave functions denoting adiabaticity . 

The effective potential energy surfaces governing nuclear motion 
in this representation are 

uJj^(R) = Wj^(R) . (Ib.l4) 

The only nonzero couplings in (Ib.9) are the nonadiabatic inter- 
actions T^ j and T^ j , hence the designation of adiabatic represen- 
tation. Figure 1 shows a perspective view of the electronically 
adiabatic surfaces of the H'*'D 2 system with a translational degree 
(R) and a vibrational degree (r) of freedom [R=(R,r)]. 





Hence in matrix notation 


T' 


= I 




k \ 



(Ib.l7) 


Similarly T" can be written as 


T" = I (P 

■ k - 


(k) si-(k) 


+ p 


(k) 




(Ib.18) 


Comparison with. (Ib.l2) shows that (Ib.18) is indeed equivalent 
to it. Denoting the colximn matrix by X/ (Ib.9) can be 
rewritten in matrix form as 





1 

2y, 


■(p 


(k) ^(k) 


+ p 


(k) 


^Oc), 


+ J_ -(k) j 

\ - 


+ D(R) - 1E]X(R) 


(T + y(R) - lE))j(R) = 0 , 


(Ib.l9) 

where 1 is the unit matrix. We see that T can be identified 
as the generalized kinetic energy operator of the system. It is 
useful to write T as 


T-= I 

k 


_l_^(k) ^(k) 


(Ib.20) 


where 




+ Ip 


(k) 


(Ib.21) 


is the generalized^momentum matrix of the k nuclear coordinate. 
Thus the matrices P and U (R) give complete dynamical infoinna- 
tion of the collision system. The diagonal elements of 
vanish identically regardless of representation. This can be 
established by applying the operator p(^) ="^ V to the ortho- 

normality condition (Ib.5) . This procedure also shows that 
p (k) is hermitian. 

Nonadiabatic couplings tend to peak near avoided crossings, 
where attains minimum values. Very often these regions 

arise because the adiabatic states are predominantly mixtures 
of two simple (molecular orbital or valence bond) structures 
(described by <J>i and 4>j) whose corresponding potential energy 
surface (Uii and Djj) cross. In one-dimensional (atom-atom) sys- 
tems, if the crossing states are of the same symmetry, then the 
"non-crossing" rule requires that the adiabatic curves formed 
from them will exhibit an avoided crossing. The generalization 
of the non-crossing rule to systems with N internal nuclear 



degrees of freedom can be stated as follows: the locus of 

points defined by the intersection of two N-dimensional exact 
adiabatic potential energy hypersurfaces corresponding to states 
of the same symmetry forms a hypersurface of at most N-2 dimen- 
sions . 

During a molecular collision process, if the relative mo- 
tion of the nuclei is sufficiently slow, they will tend to 
follow a single adiabatic potential energy surface , even near 
an avoided crossing. This is so because electronic motion al- 
ways has time to adjust to the changing nuclear configuration. 

On the other hand, if the nuclei move very rapidly in the 
vicinity of an avoided crossing, the probability of nonadiabatic 
transition (moving from one adiabatic surface to another) will 
approach unity. This is due to insufficient time for the 
electrons to adjust their motion adiabatically to rapidly chang- 
ing nuclear config^lrations . In the fast motion limit, there- 
fore, the adiabatic representation may not be the most appropri- 
ate; and it may be more suitable to choose a representation in 
which the nonadiabatic couplings are minimal . 

Representations which make use of electronic wave functions 
of desired characteristics (such as ionic or covalent) or those 
which minimize nonadiabatic couplings are in general known as 
diabatic representations - a term introduced originally by 
Lichten. Very often diabatic representations are chosen such 
that the matrices p vanish. In these representations transi- 
tions between surfaces are induced primarily by the off-diagonal 
elements Uj_j of the Hamiltonian Hgj. In the cases where these 
are small, the corresponding diabatic representation may be pre- 
ferred over the adiabatic representation, even at low collision 
energies. 

We will now illustrate the foregoing discussion with the 
simplest possible example of collision dynamics: the one- 

dimensional two-state problem. The matrix T (generalized 
kinetic energy) in (Ib.20) becomes 

T = ^ Ip^ + ^ (P.P + pP) + (Ib.22) 

where all matrices are two-by-two and P = We start with 

a diabatic representation in which P^=0 (all nonadiabatic 
couplings vanish) , and 

u'^(R) = 

To find the adiabatic representation, we have to transform 
such that 




d , . 

d . . 

'12 «> 

U22(R) 


(Ib.23) 


D^(R) = C ^(R)y'^(R)C(R) 


(Ib.24) 



becomes diagonal. The matrix that diagonalizes (R) can be 
written 


Vw* 


(R) 


cosa 

-sina 


sina 

cosa 

j 


(Ib.25) 


where 


The adiabatic 


tan 2a (R) 




'” 22 ®'-“?!®” 


potential matrix U (R) is then 


U^(R) = 


1 0 
0 1 


U(R) t 


-1 0 

0 1 


U(R) 


(Ib.26) 


(Ib.27) 


where 


tI(R) 

" h ^^11 °22 ^ ' 

(Ib.28) 

(R) 

= U^(R) + ' 

(lb. 29) 

u^(R) 

= i (U^2(R)-0 Ji(W) . 

{Ib.30) 


Of course the price we have to pay in using the adiabatic repre- 
sentation is that P^O in general. If we represent the diabatic 
basis set by the adiabatic set will be given by 


Hence 




, a “R J_ .a _ H Y ~ r 
^m£ "^m i 9R it mi 3R 


■n 


'mj 


a .,-1 

'ij 


(Ib.31) 


(lb. 32) 


e.g. 


a _ "R ^ 
^12 i 3R ■ 


We saw earlier that Pj_j_=0 in any representation and P is 
always hermitian. Therefore, P in our example can be written 



0 i] 9a (R) 
-i 0 3 r 


(Ib,33) 


P®(R) 


9a/3R can be obtained from differentiating (lb. 26) ; 

1 


a ' 


9a 

9R 


2U‘ 




(Ib.34) 


. , , a a a , a 

Recalling (Ib.22) , we require the operators P -P , pP and P p 

operating on the nuclear wave functions x (R) * We will look 

at the effects of these in turn. First, P^-P^ is diagonal. 

Therefore, it does not couple the two adiabatic electronic 

.states. Second, 



l] 

0 


’ (R) 


_ 9 _ 

9r 


(Ib.35) 


is the dominant nonadiabatic coupling matrix [which is also 
hermitian, since it involves 9/9R acting on x j ^ • Fina.lly, 



/* 

0 

-1 


l' 

0 


ki^a" (R) 


(Ib-36) 


is anti-hermitian. We will now consider the physical significance 
of the quantities a' and a" by examining the case of the Landau- 
2ener model , This model assumes straight-line potentials 
(R) and U^2 constant U^2 near the crossing region. 

From (Ib,30) we see that is also constant near the crossing 
region. At the crossing point Rq itself, U^j(R) vanishes by 
definition and in general (regardless of the choice of dynamic 
models) 

uMPq) = - U^(Rq)/[ 2U^2%)^ • (Ib.37) 


Hence in the Landau-Zener model a ’ is constant throughout the 
vicinity of the crossing point, and the quantity ARg defined by 

ARq = 1/a' (Ib.38) 

is a length characterizing the effective width of the crossing 
region, a" (R) , which governs the anti-hermitian nonadiabatic 
coupling (Ib.36), can be evaluated as [using (Ib.34)j: 


a = 


2U 


Wa"’i2>” - 


- 


u 






(Ib.39) 



which gives 


a" (R q) 




■i'(y 

W?2<V1^ 


(Ib.40) 


In the Landau-Zener model, this leads to a”(Rg)=0. Thus we can 
just concentrate on the henaitian nonadiabatic coupling term. 
Since the off-diagonal elements U ^2 diagonalized away in 
the adiabatic representation at the expense of the introduction 
of P^p, the advantage of the adiabatic vs diabatic representa- 
tions depends on the relative magnitudes of U ^2 ^nd P^p. We 
thus introduce a dimensionless constant reflecting this ratio; 


-Ea' (R)v (R) 

A.(R) = ^ , (Ib.41) 

* 2U^2<R) 

where v (R) is a state-dependent nuclear velocity given by 

X 

|vLvJ(R) = E - U^(R) , (Ib.42) 


and E is the total energy (kinetic + potential) of the collision 
system. At Rq (when v becomes state independent since 
we then have 




-Ev(Rq) 


v(Rq) 




(Ib.43) 


where 


V 

0 


-K 


(Ib.44) 


is a characteristic velocity which also plays a prominent role 
in Landau-Zener transition probabilities (see Section Ila) . Vg 
is then the critical velocity determining the suitability of 
adiabatic vs diabatic representations. If v(R)»Vg, X»l, and 
(1.41) implies that U ^2 is insignificant compared to the strength 
of nonadiabatic. couplings. Hence for high nuclear velocities, 
the diabatic representation would be more favorable. On the 
other hand, if v«Vg, u^2 exceeds the nonadiabatic couplings 
even at the crossing point, where they tend to be the strongest, 
and the diabatic representation would lose its advantage com- 
pared to the adiabatic. Far from the crossing region (R->«) , we 
can assume that 



U (co) = const, 
a 

and U*(“) = 0 , (Ib.45) 

a 

such that the logarithmic derivative of _ tends to a finite 
limit 

[U^ ] ’ 

lim ‘■12^ , ,,, 

— ^ (Ib.46) 

d 12 

^12 


Then, from (lb. 34) 


lim 

R-^ 


a ' 



lim 

R-^ 


' 

"Vl2 


12 


2U^ (“) 


and, using '(Ib.41) , 


A («) 


^v(°°) 

4U^{») 


(Ib.47) 


(Ib.48) 


Since A(“)>>1 for high velocities, (Ib.41) again shows that under 
these conditions, the diabatic representation is more favorable 
for large R. In view of the discussion above we can conclude 
that for large velocities, the diabatic representation is pre- 
ferable for all R. 


Ic. Quantum Mechanical Treatment of Nuclear Motion - Close 
Coupling 

In the last subsection we have separated electronic 
from nuclear motions by appropriately choosing Born-Oppenheimer 
electronic basis functions . To deal with nuclear motion we must 
further separate the coordinates describing relative nuclear 
motion from those describing internal motion (vibration and ro- 
tation) . Quantum mechanically this is usually done by expanding 
X£ as [in analogy to (Ib.7)] 

X^(^,R') = I (Ic.l) 

a 

... -^i 

where p stands for the relative motion coordinates and R' the 
internal motion (vibrational-rotational), coordinates. is 

the nuclear motion wave function for the a internal vibrational- 
rotational state of the Z electronic state and, as written in 
(Ic.l), is taken to depend parametrically on p. Similar to the 
case of purely electronic transitions, the nuclear wave functions 
ug^Q directly give S-matrix elements for transitions between 
electronic and internal nuclear states. To solve for we 



have to substitute (Ic.l) into (Ib.9), or, equivalently, (Ib.l9) . 
First we write (assuming there is only one internal coordinate 


R') 


T 


R 


_1L.v2 

2]i p 2u' R' 


(Ic.2) 


where u and w ' are reduced masses for the relative and internal 
motions respectively. (The case of many internal degrees of 
freedom will not affect the formal development of the coupled 
equations) . Analogous to (Ib.ll) and (Ib.l2) we have 


2 2 

i3 2\1 p 2VL R' D R 

(Ic.3) 

2 2 

% - -Ij- -1^ . (lc.4) 

Multiplying (Ib.9) by snd integrating over R', making use of 
the property (orthonoimality of for the same elecrronic 
state) 

f . 

(Ic.5) 


d\'|. 5. = 6 

la la' aa 


we arrive at the coupled equations 


- — V - E 
2v p 


U. (^ = y . , (^}u. , (^) (Ic.6) 


where 




ia,3a' 


d^R'5. (U. . + t: . t TV ,)C. , 

^la 13 1] 13 3W 


^ f- 


|d\'5. 

' la 


^ j K . , 

y y ' J Dct' 


(Ic.7) 


We note that and be orthogonal if i?^ j . The 

integral 


d rC. C. , = Y- ■ , 

la 3a' ia,3a' 


(Ic.8) 


gives the overlap between internal nuclear wave functions on 
different electronic potential energy stirfaces. The program of 



close coupling quantum mechanical calculations is the solution 
of (ic.6) , for which several powerful algorithms have been de- 
veloped. In principle (Ic.6) is capable of describing rotational, 
vibrational as well as electronic transitions- However, the 
problem may become prohibitively difficult as the nimiber of open 
channels increases . 

Id- Semiclassical S-Matrix Elements and the Stationary Phase 
Approximation 

In this stibsection we will present a formulation of the 
classical limit (Kao) of the S-matrix due to Miller (A comple- 
mentary formulation has been carried out by Marcus.) This form 
for the S-matrix for molecular collisions takes advantage of the 
fact that heavy particle dynamics of the molecular collision 
system can be well treated by classical mechanics. We will first 
consider the case of electronically elastic (adiabatic) colli- 
sions in which the motion of nuclei follows a single electroni- 
cally adiabatic surface U^(R) . In Section III we will generalize 
our formulation to electronically nonadiabatic collisions based 
on a path integral formulism. 

For the description of adiabatic motion we can use the 
Hamiltonian 

p2 

H = ^ + h(R') + V(R« ,p) (Id.l) 

where the coordinates R' and p retain the same meaning as in 
the last subsection and P is the momentum operator canonically 
conjugate to p. V(R',p) is an interaction which vanishes as 
p-^, and 


h|n> = £ |n> (Id. 2) 

n 

defines the asymptotic (p->«) internal states |n> with internal 
energies Ejj. In a given electronically adiabatic collision pro- 
cess, we are interested in obtaining the transition amplitude 
between two internal states defined by h, i.e., the S-matrix 
element 


S (E) 


(Id. 3) 


where E is the total energy of the collision system. We will 
first establish a relationship between the S-matrix element and 
the Green's function (la. 20), which will then be computed in the 
classical limit to give an expression for the semiclassical 
S-matrix. To obtain G (E) , we first have to calculate Gq (E) , the 
unperturbed retarded Green ' s function corresponding to the 
Hamiltonian 



H - V 


(Id. 4) 


in the translational coordinate and internal quantiim number 
representation |p,n>. This is given by 

sin(k p ) ik P 

<P2,n2|GQ(E)|p^,n^> = ^ 

(Id. 5) 

pp (E-en^n-^/^ 

where k = (Id. 6) 

I L -1 


where 


2p (E-£n^) 


(Id. 6) 


and p<(p>) is the smaller (larger) of p and p^- Substituting 

1 ** j Tn\ -L .cr^ 


into (la, 18) we have, for n^f^n^, 


eo -00 


1 I '^1^1^ = dp^ dp|<p2n2jGQ(E) jp^n2> 

■*0 •'o 

X <p^n2]T(E) ip^n^><p^n^|Gg(E)jp^n^> . 

(Id, 7) 


Hence 


Pl,p - '‘’2"2l= »=>|PiV = p^.p^^pj * ^V2> 


X <k2n2|T(E) |k^n^> (Id. 8) 


,00 .00 , sink2p2 

where <k n jT|k n > = dp' dp' <p' 

. -*- ■*■ Jq -^Jq ^ 


<p^n2|T(E) |p^n^> 


sin k^p ■ 


(Id. 9) 


Since the T- and S-matrices are related by 


^,ni ^^n^n^ - 2i[^ (k^k2)^^^<k2n2jT(E) [k^n^> , 


(Id. 10) 


we can identify the S-matrix in terms of G (E) as (n^rn 2 ) 



S (E) = - (k k ) exp (-ik p - ik p ) 

U 12 Pi'p2'^ ^ ^ ^ ^ 


=< <P2n2|G*(E)|p^n^> = 


^2, , 1/2 


'■ V 


X exp(-ik^P3_-ik2P2) 


dt exp 




j 




''-iHt 


^ mPi^1> 


where we have made use of (la. 20) . 

So far the treatment has been exactly quantum mechanical. 
In order to evaluate S in the classical limit, we use the 

^2'^1 

classical Hamiltonian in (Id. 11) , which is given by [compare 

(Id.l)] 


H(P,p,n,^ "" 2^7 + V(p,n,^ . (Id. 12) 


This Hamiltonian has been expressed in terms of the action-angle 
variables {nj_} and {qi}/ i=l/...,N-l, where N is the total 
number of nuclear degrees of freedom. These variables describe 
internal motion classically and, in particular, the action vari- 
ables ^ are the classical counterpart of the quantum numbers for 
the inteinial degrees of freedom. WKB quantum conditions require 
that n^ be integers at asymptotic regions (p-^) , and this forms 
the basis for. the semiclassical description of the initial and 
final quantum states of the collision system. 

In order to proceed further, we have to digress to pre- 
sent a result for the classical limit (li-M)) Green’s function 
in the cpordinate representation: , 


<q^ ] exp 



Irr > - 

(q2» 

.-K 

IV - 

2ritr 


-il/2 


X ejq) 




(Id. 13) 


where <j> is the classical action, the time integral of the classi- 
cal Lagrangian , 


<j> (q. 


- 1 


dt{p(t)q(t) - H[p (t) ,q(t) ] } 


(Id. 14) 



and p,q are the canonically conjugate momentum and coordinate. 
The only approximation involved in the derivation of (Id. 13) is 
the stationary phase approximation fomula: 


r 

if (x) 


2TriK 

1/2 


dx g (x) exp 

_ _ 

- 1 .(V 


exp 

_ di 


(Id. 15) 


where the points of stationary phase are given by 

f (x^) = 0 . (Id. 16) 

Since (Id, 15) becomes exact as-dT-K), one may say that classical 
mechanics is the stationary phase approximation of quantum 
mechanics. 

Making use of (Id, 13) in (Id. 11) we have 


S (E) = 

Pl'P2"^ 








f" ;t)/3 {p n )3 (p n 

X at ^ =—= — 

•'0 ^ (~2TriS)^ 


X exp{^[Et + (}>(p 2 n 2 ,p^n^;t)]} (Id. 17) 

where the action integral in the mixed representation (p,a) is 

ft^ 

(fi (p 2 n 2 ,p^^^;t) = I dt' [P(t*)P (t') -q(t')-n(t') - H(t')] , 


(Id,18) 


" "=2 - ^1 ' 


and the momenta and P 2 r with directions accounted for, are 
given by 

?! = ; 


P2 = -'^^2 • 


The stationary phase approximation (Id. 15) is now invoked again 
to do the integral in (Id. 17} . First we note that 



3(j> -X -i- 


E (Id. 19) 


where E (p 2 n 2 ^ p ,n^; t) is the total energy for the trajectory with 
the indicated boundary conditions. It is distinct from the con- 


stant quantity appearing in the argument of ^ 
the stationary phase condition 2 1 


(E) . In fact. 


(Id. 20) 

defines the functional relationship between E and t. The sta- 
tionary phase in the integrand of (Id. 17) can now be written 


■(n^/n^fE) = - P 2 P 2 + Et - dt’ [P (t ' ) P (t ' ) - q(t’)-n(t') 


-H(t')] = dt[-p(t)P(t) - q(t)*n(t)] , 


(Id. 21) 


3E 3(P2ff2)3(p,n^) 


i§ (n2,n^;E) 


(Id. 22) 


The pre— exponential derivative in (Id. 22) is actually the de- 
terminant 


8ni3n2 


9n23Pi 


(Id. 23) 


3P29ni 


3P29Pi 


which can be shown to be (with some messy algebra) : 


3^<j> _ u 3t 


^ fstr" ’ *(n2,n^>E) 




(Id. 24) 



The expression (Id. 22) for the semiclassical S-matrix then 
simplifies to 


^ (E) = 

n2,ni 


9 $ {n2/n^;E)/9n29n^ 

, o •4 ,xN-1 

_ '(-2irxil) 


i$ (n2,n^;E) 


(Id. 25) 

We pan go one more step to simplify the pre-exponential factor. 
From (Id. 21) we have 

<j>0t2P2,Vi;t) = *(^2,^3_;E) + ^2^2 - Vl ■ 


(Id. 26) 


and from (Id. 18), 


9n, % 


(Id. 27) 


Differentiating (Id. 26) with respect to n^, we obtain 


9$ , .V 1 

(n2,n^;E) = % + P^_ 


= ^1 "^f’l 9^ 


Up^ 9e(n^) 


^1 P 


1 1 


= q’ (Id. 28) 
9n, 1 


where in the last equality we have used (Id. 6) . If we consider 
the time dependence of qj_ without the interaction potential in 
(Id.l2) , i.e., the time dependence in asymptotic regions, we see 
that (using Hamilton's equations) 


"^1 9n, 9n, 


(Id. 29) 


■ _ 9H _ ^1 

‘’l 9P U 


(Id. 30) 


Hence 


q(t^) = const + 1^ 


(Id. 31) 



(Id. 31) ■ 




n(t 2 ) = ^2 ~ specified integers' 
p (t^) = asyn^totically large (Id. 35) 

P(t2) = {2p[E - e(n2)]}^^^ 
q(t 2 ) = anything 

If there is more than one trajectory fitting the prescribed 
bondary conditions, (Id. 33) must consist of a sum of terms, one 
for each such trajectory. To evaluate the derivatives Sn^/Sq^' , 
we must integrate Hamilton's equations to obtain a function 



n 2 (q£rni;E), whose values, of course, need not be integers. We 
will then have to solve 


^2 (Id. 36) 

(with the arguments and E of n 2 suppressed) for the initial 
angle variables which would lead to the prescribed final 
integral values for ^2- (Id. 36) is a set of N-1 equations in 
the N-1 unknowns q£, and the pre- exponential factor 9'n2/3qi is 
in fact a determinant of partial derivatives, which are to be 
evaluated at the roots of (Id. 36) . Each set of roots will 
then correspond to a distinct trajectory and will lead to a 
different ^ . 


II. SPECIFIC THEORIES AND DYNAMICAL MODELS 


In this section v?e will begin by studying simple dynam- 
ical models for nonadiabatic electronic transitions. We will 
consider low energy collisions such that the characteristic ex- 
tent of nonadiabatic coupling is small. Under such circum- 
stances the potential surfaces and the nonadiabatic coup- 

ling matrix elements can very often by approximated by analytic 
functions of R. Moreover, the transition can be considered as 
due solely to the coupling between two electronic states , thus 
the general multistate problem defined by (Ib.9) can be reduced 
to the relatively simple two-state problem. We will state our 
problem in terms of the semiclassical time-dependent formulation. 

The Schrodinger equation to be solved is 

Hii;(R,r;t) = ill ~ (R,r; t) (II. 1) 

ot 


where [compare with (Ib.7)] 


If {<|)£} 
comes 


where 


i;;(R,r;t) = I (R; t) <j)^ (R,r) . (II. 2) 

is chosen to be an adiabatic representation, (II. 1) 




3r' 


j?^i 
2 


(II. 3) 


T. . = Of*. |- W ^ •> - — <(j) I-^ !!).> — , 

13 i' 2y g^2 ^3 p ^i'3R ^3 3R 


(II. 4) 


be- 


and a one-dimensional problem has been assumed. [It is seen that 
= T|j + T^j, on comparison with (Ib.ll) and (Ib.l2) ] . The 



Tj_j are/ of course, the nonadiabatic couplings discussed in 
Section Ib.2. 

So far, our time-dependent formulation is rigorously 
quantum mechanical, and {II. 3) is equivalent to (Ib-9) . The 
semiclassical viewpoint is introduced in assuming a classical 
trajectory 


R = R(t) - 


The electronic Hamilton then becomes time-dependent through 
R(t), and the electronic motion can be described by a time-de- 
pendent wave function (r"; t) satisfying 


H , (r,R(t))$(r;t) = ±K~$(r;t) . (II. 5) 
ex ot 

* can be expanded in terms of the ^ : 

^ ■ ft 

$(r;t) = I a^(t)ij)^(R(t) ,r’)exp{-^ U^^(R(t))dt} . 

i •' 


(II. 6) 


The nonadiabatic transition amplitudes, or S-matrix elements, 
are then determined by jajj(t)|^. Substitution of (II. 6) into 
(II. 5) gives 


3r^i J 


dt(U^, - U^.)}a. (II. 7) 
33 3 


where the nonadiabatic coupling matrix elements are 


C. . s <d> . 1-1153* ./3t> . 
ID i' 3 


(II. 8) 


For a two- state problem we have 

rt 


±Ra, 


'12 


expi| 


dt (U, 


11 


}a^ 

22 2 


(II. 9a) 


iha. 

2 




dt (U 


22 




(II. 9b) 


The S-matrix element S 2 .f.i is given by 

®2^1 " k2(t-H-“) 1^ (II. 10) 

with the initial conditions 


(-«) = 0 , 


a^ (-«) = 1 


(II. 11) 



and the coupling C ]_2 is assumed to be maximum at t=0. 


Ila. The Landau-Zener Model 


It is often more convenient to formulate collision 
problems in the diabatic rather than the adiabatic representa- 
tion (see discussion in Section Ib.2). The Landau-Zener model 
of nonadiabatic coupling between electronic states of the same 
symmetry, for instance, is formulated in terms of the diabatic 
representation of the electronic Hamiltonian, defined by 



0(Rq) - Pj^(E-R^) 


= -P^(E-R0> 

d 

= a = const. , 


where 


U(Rq) 5 U^ 3 _(Rq) = ^ 22 (^ 0 ^ ' (Ila. 2) 


Rq being the point of intersection of the diabatic curves, and 
F^/ F 2 are constants. The equation of the trajectory in the 
vicinity of Rq is approximated by a linear function of time; 

R - Rq = v^t . (Ila. 3) 

From (Ib.27) , is then given by 

= U(Rq) - j(Ri+R 2 ^ (Ila. 4a) 

U 22 = U(Rq)' - (Ila. 4b) 

where 4U (R) ={AF(R-Rq)} + 4a , AF^jF^-P^l - 

(Ila. 5) 


Analogous to the derivation of (Ib.32) , we obtain 

= iii9a/3t (Ila. 6) 

where a is defined in (Ib.26). Using (Ib.34) in which ' means 
3/9t) , (Ila. 6) reduces to 



(IIa.7) 


ifia (AF) Vq 

Q _ 2 ^ ■ 

• {APCv^t)} +4a^^ 

The exact solution of (II. 9) using (IIa.4) and (IIa.7) has been 
carried out by Zener, and it is found that, for a 'once passage' 
through the nonadiabatic coupling region, the S-matrix is given 

by ■ 

As discussed in the previous section, it is thus seen that when 
the velocity Vg is small, the nonadiabatic transition probability 
is also small. On the other hand, for vg large, 32-^1 tends to 
unity, meaning that the colliding particles tend to follow the 
diabatic path. 

Ilb. The Demkov Model 

In this model (R) is given by 

0^2 (R) const , ~ ' (11b. 1) 

U^2 ^ , (lib. 2) 

and the straight-line trajectory R(t)==vgt is also employed. 
Transforming to the adiabatic representation, we have 

U^l(R) = - |[A^ + 4A^exp(-2XR) (lib. 3a) 

U 22 (R) = + 4A^exp (-2XR) ] . (Ilb. 3b) 

At asymptotic regions (R) (R) and (R) u| (R) . 

Hence in this case, opposite to tAs Landau-Zener model, the non- 
adiatic transition amplitudes are the same as the diabatic ones. 
Expanding $(r;t) in the diabatic basis set {<{>^}: 

* ' t 

®(r‘;t) = ^ b. (t)ij>f (R(t) ,r‘)exp{- dtuf (R(t))dt} , 

1 

(Ilb. 4) 

we obtain, analogous to (II. 9), the set of coupled equations 
H^i - , (IIb.5«) 




lie. Stuckelberg Theory - Curve-Crossing Problems 


In the diabatic representation (in which all nonadi- 
abatic couplings vanish) the coupled equations of motion for the 
one-dimensional two-state problem are 



^1 ^ 

(lie. la) 


+ k^(R)X2 = a(R)Xj_ 

(lie. lb) 

where 

^,2 = 3 ' 

-n 

(lie. 2) 


2p „d 
a = — ^ D 

2 12 

(lie. 3) 




and "denotes second derivative with respect to R. (IIc.l) is 
often used to describe collision dynamics with curve crossing, 
i.e., when (R) and U^ 2 intersect at some point. The 
Stuckelberg method to solve these equations is based on semi- 
classical phase integrals and their analytic continuation in the 



complex R plane. For unit amplitude entrance on channel 1 
(IIc.l) must be solved siibject to the boundary conditions 


ik (“)R ilc (“)R 

® - ^11 


^2 


- ^12 " 


ik^ (“) R 


(lie. 4a) 
(lie. 4b) 


where S]_j_ and Sj _2 are the scattering matrix elements. 

The. procediire of Stuckelberg to solve the equations is 
as follows. First one eliminates Xo (IIc.l) to obtain a 

foiarth order differential equation for xi * ^ solution of the 

form 

Xt = exp{^ S + S +"tiS + } (lie. 5) 

1 Ta 0 1 2 


is then substituted -into the fourth order equation to obtain 
equations for Sq/Sj^/.... when coefficients of the various powers 
of-^ are separately equated to zero. Keeping only Sq and Sj_ one 
can arrive at the general WKB solutions for Xi and X 2 * ^^e 

solutions thus found are observed to have singularities at R* 
which satisfies 





0 . 


(lie. 6) 


R* is in general complex and (lie. 6) implies that R* is also a 
solution. The WKB forms X]^ and X 2 then cannot be used in the 
neighborhood of R*, including the immediately adjacent real axis. 
It is, however, desirable to consider the approximate solutions 
Xi 2 analytic in sufficiently remote regions from the con- 
jugate pair R* and R* so that Xi 2 R»ReR* can be joined 
smoothly to Xi 2 R 3 _,R 2 <<R<<ReR* (Ri,2 being turning points) 

by tracing the WKB solutions along a contour remote from the 
neighborhood of R* and R*. Any contour chosen will Invariably 
pass through the so-called Stokes lines, which are defined by 


(•R 



Re 

dR(k^ - =0 

(lie. 7) 

•'r* 



1,2 

-K 

(lie. 8) 


• ^ 
and Wj^ 2 '*^be adiabatic curves obtained by diagonalizing U . 

The WKB solutions consist of dominant and subdominant solutions 

in the complex R plane. As one crosses a Stokes line the co- 



efficient of the subdorainant solution must change discontinuous ly 
in order for the solutions to join. This is known as the Stokes 
phenomenon, and one can say that a nonadiabatic transition is the 
physical manisfestation of the Stokes phenomenon. 

This connection procedure determines the coefficients 
of the WKB solutions whose asymptotic forms lead directly to the 
S-matrix. If there is only one region of nonadiabatic ity (in 
which the WKB solutions are not valid) , then there will be two 
crossing encounters due to the fact that this region is passed 
twice, once on the way in and once on the way out. Stuckelberg’s 
solution for the S-matrix is given in this case by 


S. . 
IrJ 


1,1/2 , . ^ ] 
(P . . ) exp (x(^ . . 
x,:i ^ ^x,] 


) + /^)^'^^exp(i<j> 

X , J X/ J 


II 


) (lie. 9) 


for i,j=l,2, where is chosen to be the lower adiabatic curve, 
the are the WKB phase integrals (i.e., classical action in- 
tegrals), the P's are transition probabilities; and I and II 
refer to the two crossing encoxinters and thus to the two differ- 
ent classical trajectories that contribute to the i->j nonadiahat- 
ic transition. The various quantities in (lie. 9) are given below; 



fR 

lim2(C-kR+ dR'k (R')) , 


(lie. 10a) 


lim 2(C - k.R + j dR'k (R') + 
R->« 


(•ReR* 


R 

dR'k^(R')) , 

(lie. 10b) 


♦1,2" = t,/ - - V - ’'2-' ♦ 

R-X» 


rR 

dR'k (R’) + 
^eR* 


R 

dR'k (R') 
ReR* 


j-ReR* 

+ 2 ! dR'k 2 (R')) , (lie. 10c) 


R-^ 


dR'k (R') + 
'ReR* 


dR’k^ (R') 
ReR* 


+ 2 


ReR* 

dR'k 


(R')) 


(Ilc.lOd) 



lim 2(5- k R + 

R->oo ^ 


R 


R, 


dR'k2(R') ) 


I 


(Ilc.lOe) 



II 


(Ilc.lOf) 


^2 2 = lim 2(? - k^R + 


ReR* 




dR'k^(R') + 


dR'k (R') , 
'ReR* 


11 2 
\,1 = "2,2 = ' 


^ II „ II _ 2 

"l,l " 2,2 ^ 


(lie. 11a) 
(lie. 11b) 


T, 1 T, I r, ,1 X 

"1,2 "1,2 "2,1 "2,1 ' 


(lie. lie) 


where p = ej^(-2(5) (lie. 12) 

is the local nonadiabatic transition probability and 


4: 


6 = i dR[k2(R) - k^(R)] = 2i 
Jr* 


rR* 

dR[(k (R) - k (R)] . 
^eR* 


(lie. 13) 


In (lie. 10), k£ = ki(“) and the phase factor C is given by 
(not originally given by Stuckelberg in this foiia) 


? = 65,n6 - 6 — argF (i6) - 


(He. 14) 


For 5«1, and for 6>>1, g~-l/126. 

Interference between the two terms in (lie. 9) leads to 
oscillatory behavior of as a function of E. This phe- 

nomenon is known as Stuckelberg oscillations. For example, con- 
sidering the l-*2 transition, we have 

|S 1^ = 2p(l-p) [1 + cos(?+A-B)] , (lie. 15) 
1/2 


where 


A = 


ReR* 

dR'k (R') , 

'^2 


(■ReR* 

B = dR'k^(R') . 

•^R, 


(lie. 16a) 


(lie. 16b) 


lid. Tul'ly-P'reston Surface-Hopping Trajectory (SHT) Theory 

The SHT approach is basically a classical trajectory 
approach in which nuclear motion is treated completely classi- 
cally. Classical trajectories are propagated on the entrance 



(initial) potential surface until an avoided surface crossing 
or other point of large nonadiabatic interaction is reached. A 
relation defining these avoided crossing seams have to be found 
for each problem: 


S(R) =0 . (Ild.l) 

-1. -i. 

A probability of switching to another surface P (R,R) is then 
computed and the tra j ectoiry branches . One branch , with the 
weight (1-P.) , is resumed on the old surface, whereas a new 
branch, with the weight P, is resumed on the new potential sur- 
face, with a slight velocity correction. More branching will 
occur if the branched trajectories subsequently encounter the 
avoided crossing seam again, and at each encounter, a surface- 
hopping probability has to be determined. 

The .switching probabilities at each hopping encounter 
are calculated by integrating the coupled equations (II. 9) . For 
instance, if a trajectory is initially propagated on the adiabat- 
ic surface. 1, then P = [a 2 |^ is obtained by integrating (II. 9) 
over the entire nonadiabatic coupling region. The solution ob- 
tained would, of course, depend on the choice of the classical 
trajectory R(t) used in the R-dependent quantities in (II. 9), 
which in turn depends on the choice of the potential (either 
or used for its determination. However, it has been 
found that for many triatomic systems of interest, especially 
the H 3 '*' system, the choice of the effective potential for the 
determination of R’(r) in (II. 9) has relatively little effect on 
its solution, and very often the solutions given by the differ- 
ent choices can be approximated by the Landau-Zener formula with 
appropriate parameters. Furthermore the transition probability 
depends almost entirely on the component of velocity perpendicu- 
lar to the avoided crossing seam. 

An important requirement of the model is a prescription 
for changing velocities to conserve angular momentiam and energy 
when a surface hop occurs. Since the perpendicular (to the • 
seam) component of the velocity is most effective in nonadiabatic 
coupling, velocity corrections are usually applied to this com- 
ponent only. These corrections reflect the energy defect of the 
adiabatic surfaces on the seam at the site of hop, 

Tully and Preston have obtained classical branching 
trajectories for the linear H"*" + D 2 system using the surface 
hopping model as described and compared them with trajectories 
obtained from an effective potential given by 

”eff ^ f ^1 ^ ^'^l ^ ^2 1 ^^2 ' ■ did. 2) 

If the entrance channel is on , they observed that | a 2 P , cal- 
culated by (II. 9) , was essentially zero except in the immediate 
neighborhood of the avoided crossing seam. In other regions, 
the trajectories determined from ^Qff were nearly superimposable 



on those determined from (1^) . Near the seam regions, however, 
the branched trajectories were significantly different from 
those detemined from Hence they concluded that, below a 

certain energy threshold, the probability of nonadiabatic transi- 
tions and the magnitude of nonadiabatic corrections to the poten- 
tial [given by (lld.2)] were both negligible except in regions 
near the avoided crossing seam. 

There is one important difference between the SHT model 
and the semiclassical S-matrix theory introduced earlier (and 
which will be generalized- to include nonadiabatic transitions in 
Section III) . The SHT theoiry adds probabilities (rather than 
probability amplitudes) for different trajectories and hence is 
incapable of producing interference effects leading to resonances. 
This defect, however, may be unimportant when a large number of 
initial vibrational arid rotational phases have to be averaged 
over, since the averaging process would serve to wash out the 
interference effects. 


lie. The Meyer-Miller Treatment of Classical Electronic Degrees 
of Freedom 

Even though semiclassical methods have been largely 
successful in the treatment of molecular collision processes, 
their formulations are plagued by a fundamental shortcoming - 
electronic and nuclear degrees of freedom are not treated on the 
same dynamical footing. Since it is never seriously doubted that 
classical mechanics can be applied safely to nuclear motion, 
some attempts have been made to redress the situation by advanc- 
ing schemes to treat electronic motion by classical degrees of 
freedom also. These attempts by no means aim at reducing elec- 
tronic motion to the tempting framework of classical particles 
orbiting around the nucleus. One has rather to search for the 
right classical dynamical variables which will describe the 
relevant aspects of electronic motion in molecular collision 
dynamics . For this purpose we will present a recent formulation 
due to Meyer and Miller. 

To motivate the choice of a suitable set of classical 
electronic degrees of freedom we will f irst-formulate the general 
N electronic- state problem in the diabatic representation. The 
semiclassical description introduced at the beginning of Section 
II [assigning a specific classical trajectory R=R(t) ] leads to 
the set coupled equations 

“i " ”4 ^ a. (t)U^. (R(t)) , i,j=l, ,N (Ile.l) 

for the expansion coefficients of the electronic wave function 



(lie. 2) 


®(r,R(t)) = ^ a^(t)4>^(r,R(t) ) 
i 

which satisfies the time-dependent electronic Schrodinger equa- 
tion 


H T (r,R)?i = iKB'i’/St. (lie. 3) 

el 

[Compare (IIb.5) . Knowledge of cti(t) then gives directly transi- 
tion amplitudes for electronic transitions between the states 
(In this s;:±)section we limit ourselves to one-dimensional nuclear 
motion . ) 

One of the ways to assign the trajectory R(t) in (Ile.l) 
is to assume that nuclei motion follows a force field determined 
by Ehrenfest's theorem: 


8H (r,R) 

p E pR(t) = - <*| — |<J> 

= - (t)d, (t) ~ (R) . (lie. 4) 

Since the a's are unlcnowns in (He. 4) , this equation has to be ad- 
joined to (Ile.l) to form a set of N+1 equations for the solution 
of the N+1 unknowns ai(t) and R(t) . Other choices of trajectories 
can be made (for example, as in the SHT model, in which R(t) is 
determined by the adiabatic potentials or Wgff) , but in any case, 
a single unique classical trajectory is used which is independent 
of the final electronic state, a clearly dynamically inconsistent 
situation. (In fact, this is one of the formal weaknesses of the 
SHT approach) . 

The Meyer-Miller approach consists in searching for a 
set of classical electronic degrees of freedom which would lead 
to Hamilton ' s equations, of motion completely equivalent to 
(He. 2) plus the Ehrenfest trajectory. Yet the boundary condi- 
tions of these Hamilton's equations, unlike those of the time- 
dependent semiclassical treatment, would be dynamically consis- 
tent, since electronic transitions now correspond to different 
initial and final values for certain of the classical electronic 
degrees of freedom, and different transitions entail different 
classical trajectories. These transitions can be conveniently 
and unambiguously described by the semiclassical S-matrix formu- 
lation of Section Id. 

The sought-for classical electronic degrees of freedom 
are related to the a (t) by 

a^(t) = / n^ (t) exp[-iq^(t)] . (He. 5) 

The nj_ and q^_ are analogous to those introduced in Section Id for 



the description of internal nuclear motion and can be correspond- 
ingly considered as electronic action and angle variables. The 
classical Hamiltonian Hg 3 _(n,q;t) is then given by 

Hgi(n,q;t) = <0 (r,R{t))j 


= y a*a .U*^ . (R (t) ) 
^ 1 1 i: 


i: 


(II. 6) 


that is. 


H^l(n,q;t) = I /n^.n_, exp{i (q^. -q^) (R(t) ) . (He. 7) 


ID 


1 D 


X ' ij 


Hamilton ' s equations for n^ and q^ are then given by 


% = 


(n,q;t) 


3n. 

1 


(He. 8a) 


n . = - 
1 


(n,q;t) 

__ 


The full classical Hamiltonian becomes 

2 

H(p,R,n,q) =^+ y exp{i(q^-q^) (R) , 

(He. 9) 

and Hamilton's equations for the canonical variables R,p,n",^ are 


~ 9P y ' 


(He. 10a) 


^(t) = - 11= - I ^;i^exp{i(q.-q }}^uj (R) 


1 j 1 


(He. 10b) 
(He. 10c) 


n. (t) =~ = y 2/n.n. Im[exp{i(q.-q.) }u. .] . (He.lOd) 

1 9q. t in >• Jr- 

X 3 

(He. 10c) and (He.lOd) can be shown to be equivalent to (He.l) 
if (He. 5) is used in the latter. Furthermore (He. 10b) is none 
other than (He. 4) , the Ehrenfest trajectory equation. Hence the 
goal of the Meyer-Miller program is achieved by the introduction 



of the classical electronic degrees pf freedom specified by 
(lie. 5) . 


As an example we consider the two-electronic-state prob- 
lem. The classical electronic Hamiltonian (lie -7) is given by 


H^l (n^.n2,q^,q 


2 ) = ^Al ^ ^ 2^22 


Ui2Cos(q2-qi) 


(lie. 11 ) 


Making a canonical transformation to the new set of variables 


n| = n^ + n^ = 1 (conservation of electronic probability) 

= n^ , (lie. 12) 



the two-state Hamiltonian reduces to one with only one degree of 
freedom (n^/qp = (n,q) , 

H^^(n,q) = (l-n)U^^ + ^'^22 2/n(l-n) q . (lie. 13) 

However, to make the boundary conditions amenable to classical 
treatment so that n(t 2 ) = n 2 is a function of q(t]_) = q^_, we 
have to introduce the Langer-modified Hamiltonian 

He^(n,q) = (l-n)O^^ + nU 22 + 2y4 + ^ y — - n '1 • 

(He. 14) 


Electronic transitions are then represented by the boundary con- 
ditions n(tj^) 5 nj_ = 0, n 2 = 1. T 9 compute the S-matrix given by 
(Id. 33) , we have to solve the following equation for 

U 2 (q^) = 1 (He. 15) 

and each root will give rise to a separate term in (Id.33) . 


III. THE MIIiBR-GEORGE THEORY OF ELECTRONIC TRANSITIONS IN LOW 
ENERGY MOLECULAR COLLISIONS 

In Section I we have introduced the semiclassical S- 
matrix for electronically adiabatic collisions in which elec- 
tronic coordinates do not enter into the picture at all. Transi- 
tions between internal states of nuclear motion are described 
within the framework of classical degrees of freedom (action and 
angle variables) whose boundary values are quantized. In Section 



II various ways to deal with nonadiabatic transitions are con- 
sidered, some based on specific dynamical models, such as the 
Landau-Zener model, while others, such as the SHT treatment, are 
built on somewhat dynamically ill-defined boundary conditions, 
despite successes in many applications. The Meyers-Miller for- 
malism attempts to put everything on firm grounds with the bold 
step of introducing classical electronic coordinates to comple- 
ment the classical nuclear ones. This approach, though intellec- 
tually challenging and showing great promise, is as yet rela- 
tively undeveloped in its applicational aspects. In this section 
we retreat one step from the Meyers-Miller formalism, as it were, 
and discuss a treatment of nonadiabatic transitions due to 
Miller and George, which makes use of quantima mechanical elec- 
tronic and classical nuclear degrees of freedom. Though suffer- 
ing from the undesirable feature of having to mix classical and 
quantum mechanics, this formalism leads to a quite aesthetically 
appealing interpretation: all dynamics, including nonadiabatic 

transitions, is essentially described by classical motion of the 
nuclei on electronically adiabatic potential energy surfaces. 
Nonadiabatic motion is described by the analytic continuation of 
classical mechanics into complex coordinate regions. Hence all 
dynamical information is contained in the analytic structure of 
the adiabatic potential surfaces, and nonadiabatic coupling 
matrix elements, while centrally important in quantum mechanical 
treatments, do not enter explicitly into this formulation. 
Transitions between internal nuclear states are handled in the 
same way as described in Section Id. 

Ilia- The Feynman Propagator Approach 

In our formulation of the S-matrix in Section Id, 

(Id. 13 ) was the starting point. There we noted that this equa- 
tion results from the stationary phase approximation of the 
Green's function propagator <q2 [ exp [-iH (t2-ti) /E]/q^>, which is 
equivalent to the choice of a particular trajectory, q(t) , for 
the nuclear coordinate q- This particular trajectory is the one 
obtained from the solution of Hamilton ' s equations , and is thus 
called a classical trajectory. In the treatment of a problem 
involving electronic coordinates, we will be interested in the 
more generalized propagator in coordinate representation 

<q2fX2 |exp[-iH(t2-ti)/^] (Illa.l) 

where X2,x^ are the collective electronic coordinates at t2 and tp 
respectively, q denotes the collective nuclear coordinates. It 
was Pechukas who first developed the idea of building a formal- 
ism to treat molecular collisions involving electronic transitions 
based on the propagator (Illa.l) - Instead of starting midway to' 
write down the stationary phase approximation for it, we begin 
with Feynman's path integral expression: 



<q2X2lexp 





r ^2 

. r 

2 ,1 

Dq 

Dx exp-t 

‘'q, •' 

XX 

X, 


rt 


^ dt[|- \ ” V(x,q)]} , 


(Ilia. 2 ) 


where the path integrals are over all electronic and nuclear 
trajectories x(t} and q(t) that connect (xiLq^) at t]_ and (x2/q2) 
at t2, and is the electronic kinetic energy. One can imagine 
doing the electronic path integral first {although in practice 
its actual evaluation is never carried out) and write 


< q 2 ^ 2 ' 


(^] 




q x > = 
^11 


2 i 

Dq KLx^fX^jqCt) ]exp[^j 


rt. 




(q(t))''] 

where K is the electronic propagator 


K[X2,x^;q(t) ] = 


'Dx exp [exp 


(Ilia. 3 ) 


dt {T^ - V(x,q(t))}] , 


(Ilia. 4 ) 


which is a functional of the nuclear path, i.e., it depends on 
the particular choice of the nuclear trajectory q(t) . 

The times t^ and t2 are usually taken to be the infinite 
past and infinite future respectively, when q]_, q2 and the 

system is in asymptotic electronic state <{>]_{x2_) or ' which 

are eigenfunctions of the Hamiltonian Instead of 

the propagator (Illa.l) , Pechukas noted that one actually re- 
quires the propagator in 'state' rather than coordinate represen- 
tation : 


<2,q2 |exp[-iH(t2-t^)/R] 

dXi <2 j ^ 2 ^"^^ 2‘^2 I ^^ 2"^1 ^ %^ 1^^^1 | i > 


dx, 


dx^ dx^({)*(X2)<q2X2|exp[-iH(t2-t^)/fi] |q^x^><i)^(x^) (Ilia. 5 ) 


Using (Ilia, 3 ), the 'state' representation propagator takes the 
form 



<2,q_^\e.yig 


-iH 


-ti 


!l,q^> = ^Dq [q(t) ]exp[^j ^dt j ]i 

^rt -^ 4 - 


•X (q(t))^] 


(Ilia. 6) 


where 




= jdx^jdx^ij)* (x^) ^Dx exp[^| - V (q (t) ,x) }] <|)^ (x^) 

X, 


(Ilia. 7) 


is the electronic transition amplitude. It gives the transition 
amplitude (l->-2 electronic transition) with the nuclei constrained 
to follow the trajectory q(t) . S-matrix elements for the 
Inq 2 n 2 transition, where nj^ and n 2 are the initial and final 
values of the quantum numbers for the nuclear degrees of freedom, 
are then constructed in a way analogous to the development 
following (Id. 11) . 

Retracing our steps slightly, we see that the calcula- 
tion of K 22 [q (t ) ] r which is the determination of the electronic 
transition amplitude under a fixed nuclear trajectory, is none 
other than what we set out to do via the semiclassical time-de- 
pendent formalism in Section II. K 2 I none other than 


a2exp{^ 


• ft 

— dt 

I 


0^2 (q{t))} 


[see (II-6)]; and q(t) was chosen to be 


the straightline trajectory for both the Landau-Zener and the 
Demkov models . In the next subsection we will present a semi- 
classical form for K[q(t) ] which is valid for all classical 
paths. This result is then used in (Ilia. 6) , in which the path 
integral over q is done using the stationary phase approximation, 
which effectively selects out the classical paths , out of in- 
finitely many possible ones. Then everything follows exactly as 
the development stibsequent to (Id. 11) . 


mb. The Electronic Propagator and Analytic Continuation of 
Classical Mechanics 


As 

propagator s 
The K’s are 


already remarked the determination of the electronic 
K 21 (or Kj^^) is equivalent to the solution of (II. 9) • 
given in terms of the a's by: 




rt., 

1 

t. 


dt (t) ] 


K 


11 


[q(t)l = a^ (t^) exp [:i^j ^dt W^(t)] 


(Illb.la) 


(Illb.lb) 



where the W's are the adiabatic potential energy surfaces. 

We will briefly indicate the semiclassical solution for 
the two-state case as provided by Miller and George via the first 
order perturbation solution for a 2 (with .the initial conditions 
ai(ti)=l, a2(tj^)=0). (III. 9) gives 




c^^exp 


€ 


dtMW^-W^)} . (mb. 2 ) 


Hence, from (Illb.la) 




1 


■'ll 


'1 

rt 


■dt C 23 ^exp{i dt' (W -W )}exp{- 

^•'t ^•'t 

1 1 


dt'W^} 


^dt C2^exp{- i 


dt’W^(t') -~ 


'dt’W2(t')} , 


(IIIb.3) 


which can conveniently be interpreted as an integral over t of 
transition amplitudes for transitions from to W 2 at particular 
times t. One now proceeds to the classical limit by attempting 
to do the integral using the stationary phase approximation, 
i-e., one requires that 


A. 

dt 


•t 

dfW^(t') 


dt'W^ (f) 


= 0 


or 


Wi(t) - W^(t) =0 . 


(mb. 4) 


The solution of (IIIb.4) is in general complex, which means in 
(mb. 3) one needs to evaluate the nonadiabatic coupling 0^2 
a complex time. In order to avoid this, one applies the station- 
ary phase approximation to the next order of accuracy by looking 
for a time at which the phase of the integrand in (IIIb.3) is 
least rapidly varying , i . e . , 


or 



dt'W^(t') + 



■ 

dt'W2(t') 


= 0 


d 

dt 


(W^(t) - W^(t)) 



(mb. 5) 


AW = - W^ . 


where 



Expanding the phase in a Taylor series about tQf a solution to 
(IIIb.5) , one has 


■'dt'W2(t') = T - (t-tQ)Aw^ + 0 




X (AW)^ + 


(IIIb.6) 


T = j dt (t) t J dt (t) , 

^1 ^0 


AWq = AW(t=t^) , 


(AW)" = (d^(AW)/dt^)^ ^ 
u t— 


(IIlb,3) becomes 

K2i[q(t)l = - i (C^^) at exp[i4W^(t-t„) 

0 ^ 


i (AW)5(t-t^)3j (IIIb.7) 


where the integration limits have been extended to ±<». The 
integral in (IIIb.5) is recognized as the integral representa- 
tion of the Airy function, and invoking its asymptotic form 
gives 

K [q(t)]= - <(^ I— 9 > — 

1 3t 2 t^ (2AW ./(AW)")^'^^ 


exp - 


^ AW„^2AW„ 
2 0 0 


3 -E LWW)^ 


- i °dt W, (t) - i ^dt W^(t) . (IIIb.8) 

-BJ 1 -lij 2 

1 0 


This approximate solution is recognized to resemble the exact 
Stuckelberg solution (WKB phase integral solution) for (II -9) 
most closely if the pre-ejqjonential factor is set equal to 1. 
We thus write , 

ftr, ^ ft, 'I 


K 2 i[q(t)] = e exp 


f, ft . ft 

■[e| "-E ” 2 ^^^ (IIIb.9) 

t- t _ 



1/2 


where 


6' = -r 




3 


2AW, 


(AW) ” 


(mb. 10 ) 


but whose general form is as yet undetermined. This general form 
can be obtained as follows. We already remarked that (IIIb.4) 
may not have real solutions, but it will always have complex 
solutions. For the avoided curve crossing case where tg is the 
solution to (mb. 5) , we may expand AW about t^ to obtain an ap- 
proximate solution for (mb. 4) : 


AW(t) = AW^ + O + ^(t-t^) ^ (AW ) " . 


(mb. 11 ) 


Hence t^, a solution to AW=0, is given to the above approxima- 
tion by 


= t^ + i 
* 0 


2AWq 1 


1/2 


(AW) ^ 


(mb. 12) 


and it is immediately seen that t* is also a solution. One then 
obtains • . 


dt AW(t) = 


•'t* 

* 


>4-* 




= i - AWo 


(AW) y 


= 2 


dt AW(t) . 


(mb. 13) 


Comparing with (Illb.lO) we see that 

ft^ 


6 = - — 
-ti 


dt AW (t) , 


'0 


(lllb.l4) 


which is the sought for general fomn, independent of the approx- 
imation introduced for AW(t). (mb. 9) can then be written as 


K 2 i[q(t)] = exp 


- if“dt - if*dt w^(t) - 
1 0 


'^r 

■t 


dt W^(t) 



(IIIb.15) 


= exp 



dt [q(t) ] 


ift. 


dt 


[q(t) ] 




This last equation gives the general semiclassical solution for 
the electronic propagator. 

With (IIIb..l5) , the propagator for nonadiabatic transi- 
tion (Ilia. 6) becomes 


0,q2iexp 


-iH {t^-t^)'^ 




^2 i 

Dq exp[^ 


: ft 


2,^ 1 ,• ,4.^ . 2 i 

dt j p(q(t)) 


dt Wj^(q(t) ) 


- 

'^Jt, 


dt W2(q(t))]. 


(IIIb.l6) 


It is instructive to compare this equation with the expression 
for the electronically adiabatic propagator 


<q^ 1 exp 


-iHCt^-t^)- 


-K 


IV 


“Dq exp 




(q(t))' 


-ip 


dt V(q(t) ) ] , 


(mb. 17) 


where V (q) is the single potential surface for the adiabatic 
motion. (IIIb.l6) has the following physical interpretation. 
Instead of moving on the single surface V (q) , the nuclei make 
the nonadiabatic transition at the complex time t^^ from one 
adiabatic surface to another (Wj^-^W 2 ) . Since t^ is complex, non- 
adiabatic transitions are in general classically forbidden pro- 
cesses, and the interesting conclusion emerges that classically 
forbidden processes can be described by the analytic continuation 
of classical mechanics into the complex time (or coordinate) 
domain . 

As mentioned earlier, the nuclear path integral in 
(mb. 16) (which, like its electronic counterpart, is never 
evaluated in practice) is carried out again by the stationary 
phase approximation to give; 


<2q2|exp 


piH(t -t )^| 

2 1 M 1 _ ^ 

N 

f8^2] 


-1/2 


-h J I 


UpJ 


©xp 

-ir 


(mb. 18 ) 



where N is the number of nuclear degrees of freedom, p^, is the 

momentum conjugate to q , and is the action integral 

X ^ / X 


♦2,1'%'%' ■ 


f^2 1 , 2 

f^* 

Mt - }iq(t)^ - 

dt W^ - 


'dt W„ 


(mb. 19) 


This integral is evaluated along the classical trajectory (ob- 
tained by integrating Hamilton's equations) which is determined 
by double-ended boundary conditions qj_ and q _2 which change 
from potential W]_ to W 2 at complex time t*. 

Analogous to (Id, 33) , the ej^ression for the nonadiatic 
transition S~matrix element can be written 




i ( 


i (-2TriK) 
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[aql'J 

V 


- 1/2 


X exp 


-tl 


dt{p(t)P(t) +q(t)n(t)} 


} 


(IIIb.20) 


where p,P are the nuclear translational and q,n the nuclear in- 
ternal degrees of freedom. (Note that in previous equations of 
this section g stands for the totality of nuclear degrees of 
freedom) . In the evaluation of the integral in (IIIb.20) , the 
complex trajectories will yield exponential damping factors 
characteristic of classically forbidden processes. 

Finally we point out that in the case of one-dimensional 
nuclear motion (atom— atom collisions) , (IIIb.20) reduces to the 

Stuckelberg solution given by (lie. 9) to (lie. 13) , with 5 = tt/4. 

IIIc. Examples of Nonadiabatic Transition Probabilities in the 
Miller-George Theory 

We will consider the Landau-Zener and the Demkov models 
in this subsection. The Landau-Zener model is given* in the 
diabatic representation by (Ila.l) and (IIa.2) . In the adiabatic 
representation the potential surfaces Wj^ are given by (IIa.4) , 
(here stands for U^^^) . To calculate the local transition 
probability e"^^ [see (IIIb.7)], we need 6 as given by (IIIb.l2). 
This can be changed into an integral over R, the translational 
nuclear coordinate,- if we assume the straight-line trajectory 
(IIa.3) . We have, for a 'once passage' through the nonadiabatic 
coupling region: 

^2,1 " l®2-^l!^ " exp (-26) ■ (IIIc.l) 



where 


- - — f"* 


dt ,AW (t) 


2fiv 


R. 


dR AW(R) , 


(IIIc.2) 


O-'R* 


and R* = R(t=t*) is the complex coordinate at which AW=0. From 
(I3;a.4) , 


Hence 


2 22 1/2 
AW = [(AF) (R-Rq) + 4a j 

R* = Rq + 21 a/(AF) . 


(IIIc.3) 

(IIIc.4) 


To do the integral in {IIIc.2) it is convenient to change the 
variable of integration to x, where 


R = Rq + 2i X a/(AF) . 


(IIIc.2) then becomes 


6 = 


2a 


fl 


■ITvqAF 


dxv^-x^ 


(IIIc.5) 


Elementary integration yields 

2 


5 = 


ira 


"ITvqAF 


Finally 


/ -. I* \ r"*27ra. , 

^2',l = 


(IIIc.6) 

(IIIc.7) 


which is the same result as obtained in (IIa.8). 

We now consider the Demkov model given in the adiabatic 
representation by (lib. 3), which gives 


AW = [t? + 4A^exp(-2lR) 


(IIIc.8) 


The solution to AW=0 gives a sequence a complex roots 


R^ = - + |J(2n+l) , n=0,l,2, 

(IIIc.9) 

To do the integral (-IIIc.2), we use the principle root, i-e., 
the one which is closest to the real axis and leads to the 
largest F- • We pick 

2 / -L 



(IIIc.lO) 


- A " S • 

Effecting the change of variables R->x, where 

- X ^ • 

the R-integral in (IIIc,2) becomes 

rl 


ttA . 


2Anrv 


, -iirx, 1/2 
dx (1+e ) 


0 •'-1 


(IIIc.ll) 


(IIIc.12) 


This gives 


5 = 


ttA 


2A-Kv„ 


and 


= expE-irA/dhv^) ] . 


{IIIc.13) 


For large 6, this is essentially the same result as obtained in 
(IIb.7) . 


Illd. Analytic Structure of Adiabatic Potential Surfaces 


The most appealing nature of the Miller-George theory 
lies in its description of nonadiabatic transitions; a non- 
adiabatic transition can be described by a classical trajectory 
moving continuously from one adiabatic surface to another, which 
become degenerate at some complex time (coordinate) point. In 
fact the two surfaces can be pictured as two pieces of the same 
surface, so that what appears to be a transition in real time 
(or coordinate) space is really no transition at all when the 
trajectory is allowed to be complex. The dynamical information 
embedded in the nonadiabatic couplings, which are not required 
in the present formalism at all, is completely contained in the 
analytic structure of the adiabatic potential surface, i.e., 
when is regarded as a complex function of either R or t. This 
is indeed an intellectually pleasing and physically cogent way to 
look at nonadiabatic transitions. 

The theory of analytic functions provides the natural 
framework for the description of the adiabatic potential siarfaces 
For the two state case 


W, 


1,2 




2 


(Illd.l) 



The surfaces 2 roots of the quadratic equation 

(uJ^-W) (U 22 -W) - ° ' {IIId,2) 

and as such are the two branches of the same analytic function 
w (R) , with the branch point at R^ where 

AW(R^) = ° (Hid. 3) 

W-]^ and V ?2 are also known as the two Riemann sheets of the same 
analytic function W(R), which is double- valued. W can, of 
course, also be considered as a function of t. Let the complex 
variable z be either R or t. Starting from an arbitrary point 
Zq and following any contour which winds around z^ (the branch 
point) once and ends up back on Zq, or traverses a branch cut 
(emanating from zq) once, one would go from one Riemann sheet 
of W to the next. Any contour, however, which does not encircle 
any branch point or which does not traverse any branch cut will 
force W{z) to stay on the sheet it begins with. It is at the 
point where the contour traverses the branch cut that goes 
over to W 2 or vice versa. Consider the first equality in 
(IIIb.15) . The sum of integrals can be written as one complex 
integral 

dt W(t) , (IIId.4) 

where the contour, from t^ to t 2 , is represented in Figure 2 
(solid line) . Because of the analytic natture of W, the contour 
does not have tO' go exactly through t^. In fact, any distorted 
contour from t^ to t 2 traversing the cut at some point will 
suffice (such as the dashed contour in Figure 2) . 

In much the same way we can discuss nonadiabatic 
transitions in multistate problems. Analogous to (Hid. 2) the 
N adiabatic surfaces of an N-state problem are obtained by the 
solution of the secular equation 

det[u - Wl| =0 . 


(Hid. 5) 



The NxN determinant leads to an Nth~order polynomial in W. The 
N different roots of this polynomial are essentially rhe N dif- 
ferent branches/ or Reimann sheets, of the Nth root analytic 
function W. As in the 2-state case, trajectories that make 
transitions from one branch to another can be really considered 
to be propagated on a single surface, although they would be 
required to go around appropriate branch points . 

For a transition from state a to 6 in an N-state prob- 
lem, the nuclear path integral can in general be written [com- 
pare with (IIIb.18) and (IIIb.l9)]: 


<3q2jexp 
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p(q(t))^ - W(q(t))]} (IIId.6) 



Fig. 2. Contour of integration ( — ) for the electronic propagator 
in (mb. 15) . In general to~Re(t^). The displacement is accentu- 
ated here for clarity. The dotted line represents the branch cut. 
The distorted contour (dashed) is also acceptable for the calcula- 
tion of the electronic propagator. The W's represent the sheets 
used for various parts of the contours . 



where M is the member of nuclear degrees of freedom. In (IIId.6) 

W is the N-valued analytic adiabatic surface. is a contour 

in the complex time plane connecting t^^ and t 2 and winds around 
the appropriate branch points of w in such a way that the tra- 
jectory beginning on Reimann sheet a (the ct electronic state) at 
t-j_ finishes up on Reimann sheet S (the 3 electronic state) at t 2 » 
q(t) is the classical trajectory obtained by integrating Hamil- 
ton ' s equations along the complex contour 

Before leaving our discussion of the analytic structure 
of adiabatic potential surfaces, we would like to mention a 
practical problem connected with the calculation of 6 in (IIIb.l2). 
In most cases other than the Landau-Zener model of straight-line 
potentials (for instance, the Demkov model) , there will be more 
than one branch point for W(t) . For large 6, that is, small 
transition probability, it will be usually sufficient to take the 
branch point closest to the real axis and compute 6 as though 
there were only one branch point (as we drd for the Demkov model) . 
For small 6, however, this procedure may not always be sufficient. 
For this case, the formula for the local transition probability 


p = exp(-26Q) , 

where 6 q is computed using the branch point closest to the real 
axis, may -require some correction. 

The correction procedure . usually employed is to replace 
exp (-26 q) by an expression which is the exact solution of a two- 
state exponential model due to Nikitin. The model is defined by 

d Ap a 

= B exp(-gR) + cos0exp(-gR) (IIId.7a) 

U 22 “ ® exp(-f5R) “ ^ cos6exp(-gR) (IIId.7b) 

^12 " ^21 " “ I sin0exp(-&R) (IIId.7c) 

where As, 3,A,B and 6 are all constants. U can be diagonalized 
to obtain the adiabatic potentials 

W = B exp(-gR) + ^[l-2cos0exp (-6 (R-R )) + exp (-23 (R-R ) ) 

(IIId-8) 


where R is the 'center' of the nonadiabaticity region defined by 
P 

exp(-3R ) = As . (IIId.9) 

P 


In this model Ae is the asymprotic splitting of and W 2 , and 
1/3 is the characteristic length of the nonadiabatic coupling. 
If the matrix that diagonalizes is written as in (Ib.25) , 



(Ib.26) implies that 


2a = tan 


-I 


-A sin6exp(-6R) 
_Ae-A cos0exp (-3R) , 


(Illd.lO) 


Hence a 0 at R ; (Illd.ll) 

and 2a ->• 0 at R<<R . (IIId.l2) 

P 

Assuming' a straight line- trajectory AR=Vgt, the exact solution 
of (II. 9) or (IIb.5) gives the Nikitin result for p, which in- 
corporates the effect of the multiplicity of branch points.: 


p = exp[- — 6^(l-cos0)] 


sinh [ ('!r/4) 5^ (l+cos©) ] 


sinh (176^/2) 


(IIId.13) 


This is the correction formula (sometimes known as the Nikitin 
formula) used in most applications with exponential- like diabatic 
curves. 6g is calculated by using (IIIb.l2) , with t^ being the 
branch point of W (t) closest to the real axis . 0 is obtained 

from (IIId.l2) . 


IV. APPLICATIONS OF THE SEMICLASSICAL S-MATRIX 

We will discuss two classes of problems, electronically 
adiabatic and electronically nonadiabatic problems. The former 
makes use of the formalism presented in Section Id, and provides 
the perspective and motivation for its generalization to the 
latter, which relies on Section III. 

IVa. Electronically Adiabatic Collisions 

The example we will discuss is vibrational excitation ■ 
in collinear A+BC, both classically allowed and classically for- 
bidden. Since electronic degrees of freedom are not explicitly 
considered, all dynamics is treated classically. The system has 
only two degrees of freedom, one translational R, and one 
vibrational (internal) r. The completely classical Hamiltonian 
can be written as 

p2 2 

H(P,R,p,r) = — + ~ + v(r) + V(R,r) , (IVa.l) 

2\i am 

where y and m are the reduced masses for the two degrees of 
freedom, v(r) is the vibrational potential for BC, and V(R,r) 
is the interaction potential coupling translation and vibration. 

Following (Id. 12) we have to assign (n,q) , the action- 
angle variables suitable for this problem. The action variable 



n is the classical counterpart of the vibrational quantum number 
for the free BC. The WKB quantum conditions require that 


where 



(— ) 


1/2 j- 


dr [e - V (r) ] 


1/2 


'r 

< 


£ 



v(r) 


(IVa.2) 


(IVa.3) 


and r<, r> are the classical turning points. The angle variable 
q is given by 


= (2m/fi^) 
n' (£) 


rr 


dr ' [e - v(r' ) ] 


1/2 


(IVa.4) 


where n(e) is the function defined by (IVa.2). . The Hamiltonian 
(Id. 12) becomes . 

2 

H(P,R,n,q) = ^ + V(r(n,q),R) (IVa.5) 


EHq tv 


where £ (n) is the inverse function of n(e) . For a simple har- 
monic oscillator 


/ > 1 2 , ,2 

v(r) = — mu (r - r^) 

Using (IVa.2) and (IVa.4) , one obtains 


E (n) = (n + -)'Ru 


and 


r(n,q) = + 


(2n+l)dl 


urn 


1/2 


sin q . 


(IVa.6) 


(IVa.7) 

(IVa.8) 


For the general case one can expand, the potential: 

v(r). = |mu^r-rQ)^[l + a^(r-rQ)/rQ + ^3 ^ ' 

(IVa.9) 

and £(n) and r(n,q) can then be constructed as powers series in 
A = (2ituoro2) For most diatomic X«1 and it is usually suf- 

ficient to retain terms up to first or second order in A, 

For a transition between the vibrational states n 3 _ and 
n 2 , the S-matrix, following (Id.33) , is given by 



(IVa.lO) 


(IVa.ll) 

and is the one-dimensional analog of in (Id. 34) . The 
classical trajectories q(t) and R(t) are determined by the in- 
tegration of Hamilton's equations 

q = 9H (P,R,n,q)/8n , (IVa.l2a) 

R= 3H(P,R,n,q) 8P , (IVa.l2b) 

h = -9H(P,R,n,q)/9q , 

p = - 9H(P,R,n,q)/9R . 


S_x .1. (E) = i[(-2TT±R) 




9q 


- 1/2 


. $ 


] exp[i ^ (n„,n^ ) ] 
n_ ii 2 , X 


■ij “1 


where #(n 2 ,niL) = 


dt[R(t)P(t) - q(t)n{t)] , 


We must now specify the boundary conditions for this problem. 
First E, the total energy, is an input parameter, as are ny and 
n 2 ^ the initial and final vibrational quantum numbers, both 
required to be integers . The boundary conditions are then 


n^ = specified integer 
= specified integer 
P^^ = -{2ii[E-e(n^)] 

= t{2y[E-e(n2)}^'^^ 


(IVa.lSa) 

(IVa.l3b) 

(IVa.l3c) 

(IVa.lSd) 


Note that once n^ and ng are fixed, Py and P 2 will also be fixed 
by (IVa.l3c) and (IVa.l3d) . One thus starts a trajectory with 
the initial conditions 


n^ = specified integer 


= asymptotically large 


= -{2y[E-e(n^)]} 

_ URl 

% " % * -T- 


1/2 


(IVa.l4a) 

(IVa.l4b) 

(TVa.l4c) 

(IVa.l4d) 


[which is actually the same set of equations as (Id.34) ] . In 
(IVa.l4d) £ ' (ni) iaRi/Pj_ is fixed but q^^ is as yet unspecified. 
The idea is to choose representative values of q]_ in the range 
and integrate Hamilton's equations (IVa.l2) forward in 



time until R 2 [R(t=t 2 )] is again asymptotically large and 
P 2 t(P(t))t=t ^ thus generating many trajectories with many 

final values for n 2 (which may not be integral) . In this way 
a functional relationship n 2 (q]_,nj^) is established. One then 
solves the equation 




(IVa.15) 


for q^^, where n 2 on the right is the specified final vibrational 
quantum number of the problem. Associated with each root of 
(IVa.15) for q^ is an action integral $ and a pre-exponential 
derivative On 2 / 9 qx)jj^; and the S-matrix Sj^ 2 n]_ ^ terms 

like (IVa.lO) , one for each root q^^. 

Since for any solution q^ to (IVa.15) qj^+2Tr is also a 
solution, there must be an even number of roots to (IVa.15) . 
Moreover, since n 2 (qj_) is periodic with period 2ir, any two solu- 
tions qj and (f-j-j within an interval of 2 ir must be such that 
On 2 / 8 'q’x)^ and ( 9 n 2 / 9 q 3 ^)-^^ have opposite signs. We consider 
the case wSere there are tw solutions qj and qjj in the inter- 
val 0:$q^;$2ir. Let the two action integrals (IVa.ll) associated 
with these two solutions be and respectively. Then the 

S-matrix is given by 

\n " t-2'irin|n^(q^)|]"^^^exp(i #^) + [+2iriK|n^(q^^) | 

X exp(^ * ) , (IVa.16) 

n II 


where the i outside the pre-exponential factor in (IVa.lO) has ^ 
been suppressed because it does not affect the value of j ' 

Keeping proper track of the signs and the i's in (IVa.l 6 ), 
it can be rewritten as 


'2“1 


"2^1 


1/2 .iir ^ i . , 1/2 iir i 

^1 ^ ^ii ^ ^i^ 


where 


p = [ 2 -n 1 i|n' (q ) | ] 


-1 


PlI = f27mln^(q^^)|] 


-1 


(IVa.l7) 

(IVa.18) 

(IVa.l9) 


are the pre— exponential probability factors. (IVa.l7) clearly 
reveals the quantum superposition effect inherent in the semi- 
classical formulation. The transition pro ability is 


P = [s I ^ + p.^T 2(p p )^'^^sin($ ) . 

n n, ‘ n.n^ ' *^1 II I II II I 

2 1 2 1 


(IVa.20) 


It is worthwhile to note here that the oscillatory result in 
(IVa.20) is also characteristic of the nonadiabatic Stuckelberg 



type transitions discussed in the next subsection (Stuckelberg 
oscillations ) , both being the result of quantum superposition. 
The completely classical result, which adds probabilities in- 
stead of probability amplitudes , is 


,CL 


= Px + Pii 


(IVa.21) 


So far we have only described the situation when there 
are real solutions to (IVa.15) , i.e., when the transition is 
classically allowed. However, transitions may be classically 
forbidden in the sense that no real solutions exist for (IVa.15) 
even though the state n .2 is energetically accessible [E>e(n 2 )]. 
In such cases one is forced to consider the complex roots of 
(IVa.15). To see this explicitly one may expand n 2 (q^) about 
qj^ax' real value of q^ that would lead to the maximum real . 
value of U 2 , which we denote 




max ,1 _ , 2 , 


(IVa.23) 


Suppose the designated n 2 is larger than n 2 ^^^/ so that there 
are no real solutions to (IVa.15) . Complex solutions, however, 
exist, and are approximately given by 


^1 


a 

■max 


max 


l^2^W 


1/2 


(IVa.24) 


The set of initial conditions (IVa.l4) , with complex is of 
course still good for the integration of Hamilton's equations 
(IVa.l2) . All that is physically required is that nx, Pxf ^1/ 
h 2 , R 2 / ?2 real at asymptotic regions . Q 2 is not physi- 

cally observable and can acquire any real or complex value. That 
the initial variables nx/Px'% real is specified by the 

problem.. n 2 is required to be real, which automatically implies 
that P 2 is also real [(IVa.lSd)]. R 2 , however, will in general 
be complex for large real t if qx is complex, since the 
Hamiltonian (IVa.5) couples all the variables during the propa- 
gation of a trajectory. 

Hence we seem to be running into a dilemma with the 
possible complexity of R 2 . But the difficulty only arises if 
we require t 2 to be real. If t 2 can be made complex in an ap- 
propriate way, R 2 can be made real. Actually the constraint of 
real t 2 need not be applied on any physical grounds: since in 

this formulation, E has been made exact, and by the uncertainty 
principle its conjugate variable, t, can be anything, and thus 
possibly complex. 

Suppose at an asymptotic real time t 2 , n(t)=U 2 and 
P(t)=P 2 . Propagating further in time will not change either 
n(t) or P(t). However, Hamilton's equations imply that 



(IVa.25) 


R(t) = R(t ) + — (t-t_) . 

2 jj 2 

Now R(t 2 ) may in general be complex, but R(t 2 ) can be made real 
if the final time t 2 is chosen to be 

t 2 = t^ - i ^ ImRCt^) . (IVa-26) 

This explicitly shows that in order to have final dynamical 
variables all real, the final time t 2 may have to be complex 
when classically forbidden processes are considered. 

If q-j is a complex root of (IVa.15) , it is seen from 
(IVa.24) that the other root is q|. Hence, as in the classically 
allowed case, it appears that there will also be two terms in the 
S-matrix. Now both and will be complex; but it can be 
shown that if qj is the root leading to an exponentially damping 
term exp[-Im#j], the other root will lead to an exponentially 
growing term exp[+Im®j]. Thus in the S-matrix, only one term 
should be retained, the one leading to, damping: 

S = [-2iTifin;(q^)]~^^^exp(i$^/K) (IVa.27) 

n2,ni 21 I 

and the transition probability is 

p2[lm<I. h 

^ = [2mfi|n: (q^) |] exp ^ (IVa.28) 

n^rn^ z 1 n j 

which is reminiscent of the WKB result for tunnelling in one- 
dimensional systems. 

(IVa.l7) and (IVa.27) are called 'primitive' semiclassi- 
cal results because they only work well for n 2 <n 2 ”*^^ and n 2 >U 2 ™^^ 
independently, i.e., the separate solutions cannot be connected 
smoothly as n 2 traverses n 2 ®^^* In particular, when n 2 
approaches U 2 ®^^ arbitrarily closely, both (IVa.l7) and (IVa.27) 
break down since |n 2 (qj This situation recalls the 

classical rainbow effect. To remedy it, one has to introduce 
uniform semiclassical formulas which, in addition to effecting 
a smooth connection between n 2 <n 2 ^^^ and n 2 >n 2 ”'^^, will also 
approach the 'primitive' semiclassical results in regions where 
they work well. 

IVb. Electronically Nonadiabatic Collisions 

We consider here the collinear A+BC nonreacuive colli- 
sion process again, but with electronic transitions. A typical 
process of this type can be written 



where (n 2 ) is the initial (final) vibrational quantiam number 
of H 2 / and the halogen X undergoes a spin-orbit electronic 
transition. The total Hamiltonian for this system can be written 


H = 


*2 32 ^2 32 


(IVb.2) 


where R is the translational X-H 2 coordinate and r is the vibra- 
tional H 2 coordinate, p and n are the reduced masses pertaining 
to R and r respectively. 

In the previous subsection, both R and r are treated 
classically; and the nuclei are considered to be moving on a 
single electronic potential surface v(x)+V(R,r) [see (IVa.l)]. 

The transformation (p,r)^(n,q) is then made, with the action- 
angle variables (n,q) still treated classically. Since an 
electronic transition is involved here, we have to consider at 
least two electronic surfaces (we will limit ourselves to ex- 
actly two) . There are two equivalent ways to proceed using the 
Miller-George theory. One is using the electronic representa- 
tion, i.e., finding an electronic basis set (<>(R,r,x) (where x 
stands for electronic coordinates) which will generate the po- 
tential surfaces and couplings U^. = <4i(Hel [ (the subscript 

X means that brackets denote integration over electronic coordi- 
nates only) . We can assimie that the starting basis set will 

in general be diabatic . can then be diagonalized to yield 

the adiabatic surfaces required for the Miller-George theory. 
In this approach the hypersurface W^(R,r) is two-dimensional and 
one is required to know the analytic structure of a function of 
two complex variables , very often not an easy task. The second 
way is to use what is known as the vibronic representation. This 
approach makes use of (to begin with) a diabatic vibronic basis 
set 


^in ^ <j>^(x,r,R)t^(r) 


(IVb.3) 


where ifin(t) is a vibrational v/ave function of BC independent of 
the electronic configuration of A. This set generates diabatic 
surfaces and couplings given by 


U*? . 

in,]m 


= <§. 
in‘ 


2rt 


+ H 


3r 


el 


<|) . >-i. 

x,r 


(IVb.4) 


which can again be diagonalized to produce adiabatic surfaces as 
required by the Miller-George treatment. The advantage here is 
that the adiabatic ' surfaces ' W (R) here will only be functions 
of a single variable R, and their analytic structure is much 



easier to handle than those of Wj_(R,r). Hence the classical 
nuclear problem is reduced to a on e-dimens xonal one, - and the 
result of the application of the Miller-George theory (IIIb.18) 
reduces to the Stuckelberg solution given by (lie. 9) to (lie. 14) . 
In this subsection we will use the vibronic representation. 

The DIM formalism applied to the collinear model 
(IVb.l) yields a two-by-two form for the electronic Hamiltonian: 

Uii(R,r) = v(r) + D + (R-r/2) -i- A^(R+r/2) + iA/3, i=l,2 


(IVb.Sa) 

^12 " ^21 (IVb.Sb) 

where D is the depth of the Morse well v (r) , the A^ ’ s are the 
anti-Morse functions, and A is the spin-orbit coupling of X. 

can be diagonalized to yield the adiabatic electronic sur- 
faces Wj^(R,r) with the accompanying adiabatic basis set 
such that 


H ^<{tf(r,R,x) = W. (R,r) (r,R,'^ . (IVb.6) 

el 1 IX. 


We prefer to use the representation [compare with (IVb.3) ] 


- (p. (x,r,R)i[» (r) , 
xn 1 n 


(IVb.7) 


where the vibrational wavefunctions (r) are Morse functions 
satisfying -the equation 


v(r)}i{i (r) = • 

^ ^ n n n 

<72T 


(IVb.8) 


where is the asymptotic vibrational energy for the n vibra- 
tional state of H 2 - Using (IVb.4) and first integrating over 
.the electronic coordinates x, we can write 

-^2 2 

. = 6 . .<ii 1-'^ t W. (R,r) U > 

xn,]m xj ^n' 2ri ^ 2 x ' ^m r 


k . (r,R) |i|i > - — Ofi ix. ■ (r,R)-;;^U > , 

2ri ^n'^ij ' '^m r n n'Hu '3r'^mr 




(IVb.9) 


where 


2 

?..(r,R) = <(pf(x,r,R)| — r- <j)^(x,r,R) , 

13 1 3r ^ 


(IVb.lO) 



(IVb.ll) 


X,.(r,R) = «j)^(x,r,R) 4)^(x,r,R)>^ . 

ij X or j X 

These coupling matrix elements are similar to the nonadiahatic 
couplings T|j and T^j [(IB, 11) and (IB. 12)] discussed in Sec- 
tion Ib. 

We can also show that 



2 


= I X 

13 

k=l ■ 


ik^kj “^ij' 


(IVb.l2a) 


(IVb.l2b) 


The couplings x furthermore be expressed rn terms of 

U<^: 


(3uf ./8r) 

= y r r . — 


^12 li 23 W^-W^ 


(IVb.l3) 


where 


U 


r = r 

11 22 


r = -r 

12 21 


12 




ud 

12 




(IVb.l4a) 


(IVb.l4b) 


The 'diabatic' vibronic potential surfaces can then be written, 
from (IVb.9) : 


u'^ . (R) = <i|> I-"" 

in,im n' 2ri 


u > 

, 2 i 2 b ^12 '^m r 
3r 


(IVb.lSa) 


The off-diagonal matrix elements become: 

^ 9X 


^ln,2m^^^ 2 b 


^^12 3r 


12 


3r 


liji > 

‘ r 


(IVb.lSb) 


For a particular vibrational transition n->-m, we have a 
one-dimensional two-state problem, i.e., we just have to deal 
with the 2X2 matrix (IVb.15) with fixed n and m. This matrix 
can be diagonalized to yield the final adiabatic vibronic sur- 
faces W|(R) and the S-matrix is given directly in te 2 nns of the 
Stuckelberg solution '[(IIc.9) to (IIc.l4)]. Figure 3 illustrates 
the diabatic and adiabatic vibronic curves and the coupling X]_2 
for the BrH 2 system. The potential curves correlate asymptoti- 
cally to the states 11 and 20. 




R (bohrs ) 

Fig. 3. The diahatic vibronic curves (the solid lines labelled 
Id and 2d) are split by the vibrational nonadiabatic coupling 
Xi 2 (see insert) into the adiabatic vibronic curves, labelled la 
and 2a, which are represented by the dotted lines. Note that the 
labelling is reversed inside the avoided crossing. 


As mentioned before, the formula p=exp(-26) very often 
does not yield good results for P=|sp when compared with quan- 
tum mechanical calculations; and the Nikitin correction formula 
(IIId.13) is usually employed. There is, however, no cut and 
fast rule regarding the applicability of the Nikitin formula, xt 
being exact only for exponential diabatic curves and couplings. 

In the present example, this formula has been applied to Morse- 
like potentials (IVb.5) , and found to help appreciably in bring- 
ing the results into good agreement with quantum calculations . 

Another slightly disconcerting factor in the application 
of the semiclassical approach is the difficulty with the phase 
factor C in (lie. 14) . The value tt/ 4 is mot often quoted and 
used in the literature. In the present applications, since the 
Nikitin correction formula is sued, the foimi as given by (lie. 14) 




would be rendered less rigorous , end the justification for its 
use in conjunction with the Nikitin formula requires more work. 

(In any case the use of this form is not very frequently reported 
in the literature) . In our work we have varied 5 between -ir/4 
and ir/4 to compare with quantum mechanical results. Reasonable 
agreements are achieved at either extreme, but it is found that 
5~ir/10 would shift the results into phase most closely with the 
quantum results . 2 

In Figures 4 and 5 we present some results for 2 jh1 

(written n 3 ^-^m 2 ) for X=Br. In each case the semiclassical 
results are compared with exact quantum mechanical coupled-chan- 
nel calculations. The reasonable agreement is apparent, point- 
ing to the reliability of semiclassical calculations. Similar 
results are obtained for the Br-H 2 2 j_-l 2 ^nd F-H 2 0 j ^~02 transi- 
tions. It should be noted here that the distinct oscillatory be- 
havior may be due entirely to the artifact of the collinear model. 
In three-dimensional calculations where there is strong coupling 
between translation and the various internal degrees of freedom, 
and especially when averages over a large number of quantum 
states are involved, the interference effects may be quenched. 



E (hartrees) 


Fig. 4. Transition probability versus energy curves for the 
O 1 -O 2 transition in the Br + H 2 collision system. The threshold 
for this transition is 0.027 hartree. The solid curve represents 
quantum results, and the dashed and dotted curves are the serai- 
classical results obtained with the phase factor C=0 and +ii/4, 
respectively . 




Fig. 5. Transition probability for the I 1 -O 2 transition in the 
Br + H 2 collision system. The threshold for this transition is 
0.029 hartree. The solid curve represents quantum results, and 
the -dashed curve is the semiclassical result obtained using the 
phase factor 




V. THEORY OF INTERACTION WITH FIELDS - LASER-ENHANCED COLLISIONS 


The widespread interest in lasers has generated a quest 
for a deeper understanding of the effects of intense radiation 
on molecular collision processes. In this class of phenomena, 
as distinct from its atomic counterparts, dynamic considerations 
play a major role, and field- induced molecular events do not have 
to satisfy strict resonance requirements. Thus the radiation 
field may be expected to have very prominent effects on the 
collision dynamics. In particular, such molecular processes as 
collisional excitation, unimolecular dissociation, radiative 
collisions, and Penning and associative ionization are expected 
to be strongly influenced by an intense laser field. In this 
section we will generalize the Miller-George theory to include 
radiation interaction. 

As soon as a field is put in, many problems suggest 
themselves. If we are dealing with a laser, the field may 
approach being monochromatic but always not exactly so. Even a 
single mode laser has a spectral linewidth owing to broadening 
effects such as Doppler, collision and. radiative lifetime 
broadening. How many modes of the laser field must one include 
in the interaction Hamiltonian? This immediately raises the 
question of what the best representation is for the description 
of the photon field in laser-molecule interactions. There are a 
number of candidates: the Fock (photon-number) state representa- 

tion, the coherent representation, or even the classical field 
approach, i.e., regarding the field simply as a time-dependent 
driving force on the system. What part does the polarization 
of the field play in las'er- influenced molecular collisions? 
Recently there is even some controversy regarding the fundamen- 
tal validity of using the vis a vis the interaction 

under the dipole approximation. These and other questions have 
not as yet been answered definitively. In this section we will 
restrict ourselves to the simplest possible (although not en- 
tirely realistic) situation of a monochromatic field interacting 
with a two-electronic state molecular system. The dipole approx- 
imation is assumed, and the photon-number representation is used, 
as this lends itself most easily to physical interpretation. 

Also, complications deriving from the intrinsic angular momentum 
of the photon will not be considered. Even under these limita- 
tions, the radiation field can be demonstrated to have signifi- 
cant effects on the collision system. 

Va. The Electronic (Vibronic) -Field Representation 

The field-free collision system can be described by the 
electronically adiabatic potential surfaces Wj_ (^ which satisfy : 

H^l?i(R,3j) = (R) (R,^ 


(Va.l) 



where {i^^} is the adiabatic basis set, and R (x) are the nuclear 
(electronic) coordinates. One can also work with the vibronic 
surfaces W^(R) as discussed in Section IVb, where R represents 
the translational degree of freedom only. Our discussion in 
what follows can be phrased equivalently in the electronic or 
vibronic representations; for the sake of definiteness, we will 
choose the foinner. The Hamiltonian due to the presence of the 
field is given by 

='R(j}a'^a + pT(^,x) 'E’q (w) (a’^+a) (Va.2) 

= dlwa'^a + . 

The first term represents the free-field Hamiltonian and the 
.second term, H^^, the interaction Hamiltonian under the dipole 
approximation. It is understood that the photon creation and 
annihilation operators, and a respectively, are those for the 
single frequency to; and 

Eq((o) = (2TTKio/V) (Va.3) 

where V is the quantization volume of the collision system and 
e is the unit polarization vector of the external field. 

A natural basis set to be used for the matter-field 
system consists of states of the form (ji£jn>, where |n> is the Fock 
state with n photons (of frequency to) , since they describe the 
asymptotic (R-^) configuration of the matter-field system. For 
instance one may visualize that before collision, the molecular 
system is in the i'^^ electronic state, while there are nj_ photons 
floating around in the external field (in the quantization volume 
V) , This state of affairs is described by the asymptotic state 
^i|ni>. After the collision, again at asymptotic regions, the 
molecular system may end up in a different electronic state (say 
the j) whereas the number of photons in the external field may 
be changed to nf. The final asymptotic state will then be de- 
scribed by (j)j|n£>, and the S-matrix element corresponding to this 

process may be written S. . .If n^>n, , stimulated emission 
^ :n^,in^ 3 1' 

has taken place; if nj>n;i^, absorption has occurred. 

Even though is adiabatic in the field- free Hamil- 

tonian Hgj_ , the states [ n> will no longer be adiabatic in the 
Hamiltonian Hgj|_+H^ , since the interaction Hamiltonian mixes 
them. They are only eigenstates of the Hamiltonian 

H = H . + 'Ema'^a , (Va.4) 

R el 

H (j) |n> = (W . +nilw) (|) . j n> . (Va.5) 

The potential surfaces Wj^+nKm can now be considered the diabatic 
surfaces in the presence of the field. Analogous to the field- 



free case, at locations (along it) where resonance occurs 
(W^+nifioj=Wj+nj3Taj) , avoided crossings may be produced if symmetry 
conditions permit, i.e., if 

d^x(j)^]i (x,R) 7 ^ 0 . (Va-6) 

In such cases , one may generate adiabatic curves ]ust 
as one does in the field-free situation. These adiabatic curves 
(R) satisfy the equation 


where 


(H ^ (R,x) = 

el r 1 


If; (R,x) = I a. 


£ ,n 


i£n 


E^(R)1|)^(R,X) 
(R)<i)^(R,x) [n> . 


(Va.7) 


(Va.3) 


The representation is termed the electronic-field representa- 
tion while Ei is known as an electronic-field surface. The 
Ei's are, of course, eigenvalues of H=Heq+Hf, where H is ex- 
pressed in the <|>j_|n> representation. Suppressing i in (Va.8) 
(i.e., writing it as a general eigenvector) and substituting in 
(Va.7) , we obtain the infinite set of coupled equations for 


Ea. = (W„+rtfTw) a. + I d' [/n+1 a. , + i/n a. , ] / 
£n £ £n “ £i i,n+l x,n-l 

X 

(Va.9) 


where 


d;i(R) = E'O.)- 


d\^*(R,x)t(R,x)(|)^(R,x) (Va.lO) 



<nja|0> = 0 (Va.lla) 

<n'|a|n"> = 5 , , (Va.llb) 

<n'|a’^|n''> = /n"+l 5 , , . (Va.llc) 

These properties lead to the important fact that only states 

with a difference in photon number equal to one are directly 
coupled; and hence, to first order in the dipole strength d' , 
only single-photon processes are possible when the interaction 
teirni proportional to ^ is ignored (A is the field vector poten- 
tial) . Even with this simplification, (Va.9) shows that, in 
order to obtain the adiabatic (electronic-field) surfaces, an 
infinite-dimensional matrix iias to be diagonalized in general. 



Vb. The Rotating-Wave Approximation (RWA) 

For the present discussion we limit ourselves to the 
two-state molecular system with W 2 >W^. In this molecular colli- 
sion problem, the energy detunings A+ for a transition 
(j»2!n±l> 

A^(R) = ± -hm , (yb-1) 

must be examined as a function of all the nuclear degrees of 
freedom. For W 2 >Wj^, for all % hence the transition 

<!>lln> ^ (j) 2 |n+l> can be considered anti-resonant (or potential 
energy non-conserving) . A_, however, may vanish for certain 

nuclear configurations, implying that the transition <|)x|n> 

(j) 2 in-l> is resonant at those configurations. 

We now assume that dj_x=0 (which would be rigorously true 
for a homonuclear atom-atom collision system) . In this case 
the secular matrix for (Va.9) block-diagonalizes , and it can be 
written in the general foinn shown in (Vb.2) , where the coupling 
strength is 

N is the number of photons in the external field and ^Tq is the 
electric field strength. The block diagonalization reflects the 
fact that, in making the transition system must 

gain or lose an odd number of photons . The two non-rdiagonal 
blocks are completely equivalent; and depending on the initial 
conditions of the collision environment, one can ignore the 
appropriate one. For instance, if the collision system is pre- 
pared in the state and the applied radiation field has field 
strength 


Eq = /2ttNRm/V (IVb.4) 

with N odd, then we can work entirely with the upper left block; 
and with N even, the lower right block. In practice, of course, 
N is so large that <3 n~%+i~%- 1' snd the two blocks would be 
largely equivalent. 

It is in principle possible for a system prepared in 
initially to absorb n photons (n odd ><N) or emit n photons 
(Kn odd «“) , and either stay in <j>x or make a transition to ^2* 
While processes involving transitions — ^‘i’21^2^ 

may not be resonant if n 2 <n^^, they will always be anti-resonant 
if ri 2 >nx. The RWA in the collision case (similar to the atomic 
case) consists in throwing away all direct couplings between 
diabatic states of the latter kind (n 2 >nx) ; for example, 

between the states (Ji21n- 1> and <t>ijN-2>. It is expected to be 






less accurate- when the applied field strength is very large; or 
in non-resonant situations (as opposed to anti-resonant ones) 
corresponding to such that the energy detuning,, 

W 2 ~W-|_-Kaj, is not significantly different from that for the anti- 
resonant transition, W 2 -W 3 _+'Sai. In the latter instance, applica- 
tion of the' RWA would not be internally consistent, since it 
would not be justifiable to chop only couplings connecting 
states with n 2 >nj^ while keeping those with nj^>n 2 if the detunings 
are roughly equal for the two cases. 

In our example of the two-state model, the RWA then 
leads to complete block diagonalization of (Vb.2) into indepen- 
dent 2x2 matrices. For the initial conditions mentioned pre- 
viously one has to consider only the diabatic matrix 


W^+ {N-l>truj d„ 

2 N 


N 


W^+NfftoJ 


(Vb.5) 


Diagonalization of (Vb.5) will lead to the two adiabatic (elec- 
tronic-field) surfaces 


- ®r» + ^ + 441 . 

(Vb.6) 

[Compare with (Illd.l)]. In effect the RWA limits the consider- 
ation of absorption and emission to single-photon resonant pro- , 
cesses. Figure 6 illustrates the construction of the electronic- 
field surfaces . 


Vc . An Example 


To apply the Miller-George theory, Ejl, 2 
as analytic functions of As discussed before 
garded as different Riemann sheets of the single 
tion E (^ , with branch points which are- roots of 


are considered 
they can be re- 
analytic func- 
the equation 


Ei(R) - E2(^ 


' (Vc.l) 


or 


[W^ (R) -w^ (R) --Km] ^ + 4d^ (R) = O' . 


(VC. 2) 


Classical trajectories can be propagated on E(^, as described 
before, making transitions by rounding the appropriate branch 
points . 




Fig. 6 . (a) Schematic drawing of two field-free adiabatic sur- 

faces and W2 as functions of the translational coordinate R. 

Rq is the point at which the field is in resonance with the two 
surfaces . (b) Schematic drawing of the electronic-field-surfaces 

E]_ and E2- The splitting between them at Rq is caused by the 
radiative coupling d^. 


We again consider an example of collinear A+BC collision. 
The specific process is 

Br(^P^^2^ H^Cv^O) + -Kcj (v= 0 ) , 

(Vc. 3 )' 

where v designates the vibrational state of H2 . and W2 are 
taken to be vibronic po’tential surfaces as before (see section 
IVb) so that nuclear motion is reduced to one-dimensional. The 
interaction of the field with the molecular system is approxi- 
mated by the interaction with the asymptotic collision species 
alone. In this case the magnetic dipole transition matrix ele- 
ment for the halogen dominates. Hence the magnetic field of the 




laser radiation appears for I!q (formerly considere d to be the 
electric field) in all expressions. The value’ of -/ a /3 a.u. 

(atomic units) is used for the magnetic dipole transition moment 
Il^l2i' (This value actually holds for all halogens). dTm is 
chosen to be 1.001 times the Br spin-orbit splitting (0.4568 eV) , 
such that the laser cannot induce a spin-orbit transition without 
the aid of the H 2 collision partner. A field strength Eq of 
5.1x10^ V/cm(10“3a.u.) is used. The vibronic-field surfaces 
(curves in the present case) can be constructed according to 
(Vb.6) , where 2 2 assemble the schematic drawings of 

Figure 6. 

Since this is a one-dimensional problem, the -Stuckelberg 
solution (lie. 9) to (lie. 14) again applies, with the slight alter- 
ation. that in these equations are replaced by E^. A pair of 
complex intersection points between Et and E 2 is located at 
R*=4.01±0.422i bohr (1 bohr = 0.529 A), and classical trajectories 
were integrated in the complex plane which switch curves smoothly 
at the intersection point. The Nikitin correction formula 
(IIId,13) is again employed and ? is chosen to be tt/ 4. The 
semiclassical results are compared with the quantum mechanical 
as well as the field-free results in Figure 7. The high prob- 
abilities for the laser-induced' process are in direct contrast 
to small probabilities for the same process in the absence of 
the field. The field-free process 'tends to be electronically 
adiabatic and the field is the sole agent in inducing electroni- 
cally nonadiabatic transitions in the present case. 


VI. MISCELLANEOUS TOPICS IN THE THEORY OF LASER-ENHANCED COLLI- 
SIONS 


The majority of collision processes only require a 
finite number of states (usually two) for their description. 
There are, however, some in which one or more discrete states of 
a system can interact with the continuous spectrum. There are 
in .general known as bound- continuum processes. In these pro- 
cesses one is interested in the S-matrix element connecting an 
initial discrete state and a final continuum channel; and very 
often there is no straightforward way to reduce the problem to a 
finite-state one or even one involving infinitely many discrete 
states. In this section we discuss two examples of field- in- 
duced bound-continuum collision processes, spontaneous emission 
and Penning ionization, in which S— matrix elements may be cal- 
culated using discrete-problem techniques by the introduction of 
certain approximations and discretization procedures. The for- 
malisms for both problems still leave much room for rigorous 
justification and improvement, and may not be adaptable to the 
most general collision situations. However, they do provide 
convenient descriptions, both intuitively and physically, of the 
respective processes . Our discussion of the bound-continuum 
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Fig. 7. Probabilities for the collinear reaction Br( ^2/2^ 
H 2 (v= 0 ) +-noi -> Br{2p|L^2) + H 2 (v= 0 ) as functions of initial rela- 
tive translational collision energy. The value-of-ITa) is 1.001 
times the asymptotic spin-orbit splitting. Shown are the results 
from the quantum calculations using the two- (vibronic ) state 
model (solid line) and results from semiclassical calculations 
using the two- (vibronic) state model (dashed line). The quantum 
probabilities for the field-free transition (dotted line) have 
been multiplied by 100. 


problems will be based on the Miller-George semiclassical foimial- 
ism developed earlier, but we should mention that the discretiza- 
tion procedure introduced below for the ionization problem is 
also adaptable to a quantum mechanical coupled-channels treatment. 

The semiclassical formalism we have developed in the 
present lecture series for discrete-state field-induced problems 
can be considered to be exact within the limitations of finite- 
dimensionality and the RWA. The treatment is non-perturbative, 
and collisional and radiative interactions are treated on the 
same footing. Despite these advantages, its application will 
often raise severe practical problems [e.g., in the analytic 
continuation of the potential surfaces and the location of branch 
points (seams) ] . For this reason, simpler treatments often re- 
tain their appeal and value. In particular, time-dependent 




perturbation treatments may be most useful when the laser field 
is of moderate strength and considered as a classical driving 
force. To conclude the present lecture series, we will give in 
Section vie a description of such a theory applied to laser-in- 
duced collisional processes as due to Weiner, 

Via. Collision-Induced Spontaneous Emission 

We consider non-reactive processes of the type 

A + B+iito->-A + B+ (Vla.l) 

(e.g., the Br-H 2 system) where w is the frequency of the incident 
laser radiation and o)' that of the emitted radiation. To describe 
the interaction of the collision system with the incident laser, 
we can make use of the electronic (vxbronic) -field surfaces 
introduced in Section Va and illustrated in Figure 6 . 

These electronic-field surfaces [arising from the 
diagonalization of (Vb.2)] can be regarded as fo 2 naing a spectrum 
for spontaneous emission if radiative coupling with the vacuum 
field is now considered. This interaction involves in principle 
an infinity of field modes , since a continuum of emission fre- 
quencies is possible corresponding to the different nucl.ear con- 
figurations at which emission can take place. The symmetry 
allowed emission transitions between the electronic-field sur- 
faces are illustrated in Figure 8 (both dotted and solid lines) . 

If we consider the case where the system starting on W^ , asymptot- 
ically is de-excited back to asymptotically after sequential 
absorption and emission of photons, only the transition marked 
by the solid line (in Figure 8) need be considered. 

The semiclassioal picture allows us to think of the 
photon of frequency m ' as being emitted near the configuration 
Rg; such that 


-HW = E„(R ) - E, (R ) .' (Via. 2) 

2 £ 1 £ 

(e) 

Emission-electronic-field surfaces E+ describing the dynamics 
of this event can then be constructed analogously to the elec- 
tronic-field surfaces. These are illustrated in Figure 9. It 
should be noted that the dimensions of the emission avoided 
crossing are much smaller than those of the absorption one since 
the former is generated through coupling with the vacuum field. 
Referring to Figure 9 (b) , the dynamical picture then emerges 
that if the system is still propagating on the surface E^®) at 
asymptotic regions, it will have emitted a photon of frequency 
oj' and become de-excited to the ground state W^. 

The probability of emission at R^ is determined semi- 
classically by the trajectories which avoid emission before and 
which lead to emission when that configuration is reached. If we 
consider the case of discretized allowable emission configure- 




Fig. 8 . Symmetry-allowed emission transitions in the two-state 
electronic field representation. The splitting at Rq is caused 
by the absorption radiative coupling. The solid line represents 
the de-excitatiqn transition from W2 back to Wj_. 


tions (in practice, one needs to consider a continuum) , a system 
propagating on the surface E2 after absorption will encounter a 
series of configurations at which it can emit. Let one of these 
be At this configuration, the approximation is made that 

there is emission coupling only at the frequency (Rgi) “ 

El • Couplings at frequencies corresponding to more distant 
configurations can be considered as yet not 'turned on'. Those 
at frequencies corresponding to configurations prior to R^^ are 
assumed to have negligible effects on the dynamics of the system 
at Rg-^, since the system has 'survived' pre-emission loss at 
these configurations before reaching Rej, . Hence, at each possi- 
ble emission configuration, emission-electronic- field surfaces 
can be generated independently of emission couplings at other 
configurations. A schematic representation for the discrete case 





Fig. 9. (a) Schemacic drawing of the diabatic surfaces for the 

description of single-mode emission at R^.. (b) Schematic drawing 

of emission-electronic-field surfaces for single-mode emission at 
Rg.. The avoided crossing is caused by emission radiative coup- 
ling. 


is illustrated in Figure 10. 

When the emission sites become continuously distributed 
along R, we still retain the approximation of treating each 
avoided crossing as generated by a localized emission coupling. 
In this continuous limit the discrete series of branch points 
will constitute a line of branch points in the complex R-plane. 
It will be expected to be relatively close to the real axis com- 
pared to the location of the absorption branch point because of 
the relative weakness of the emission couplings. Whereas in the 
discrete case a trajectory, after rounding each emission branch 
point, may return to the real axis before going around the next, 
in the continuum_ case it will be forced to deviate from the real 
axis throughout the configurations at which emission coupling is 
active and be propagated along the line of branch points . The 





Fig. 10. Schematic drawing of emission electronic-field surfaces 
for discrete emission configurations. is the emission 

site. 


imaginary action accumulated along this line will determine the 
accumulated probability for making local transitions from 
to continuously, that is, the accumulated probability for 

not having emitted a photon while the molecular system stays on 
E2- the actual emission site the trajectory reverts back to 
■the real axis and is propagated to asymptotic regions on . 

The shape of the emission spectrum (as a function of the 
emission configuration Rg) is then found to be determined by the 
shape factor 

P (E,R ) = exp[-2x(E,R ) ] {1 - exp [-2C (E,R ).] } (Via. 3) 

o E £ £ £ 

where E is the total energy of the collision system- The action 

factor X accounts for pre-emission effects and 1-exo (-2? ) is the 

e 




local emission probability. Since the emission branch points 
tend to depart more from the real axis as increases, 1-exp 
{- 25 g.) will increase and exp(-2x) decrease as Rg- becomes large. 
Thus it is seen that pre-emission and localized actual emission 
effects compete with each other for the total emission prob- 
ability at a particular configuration. Results of model cal- 
culations indicate, however, that this total probability depends 
critically only on the local emission probability at the actual 
emission configuration. 

VIb. Laser-Induced Collisional Ionization 

Let us consider non-reactive processes of the type 

A(A*) + B t-Km ->• A + + e"(e) PI (VIb. la) 

^ Ab"^ + e"(e) AI (VIb. lb) 

where m is the laser photon frequency, e is the kinetic energy 
of the emitted electron, and A (A*) means that A may or may not 
be in an excited state. Processes (VIb. la) and (VIb. lb) are the 
counterparts to field- free (FF) Penning ionization (PI) and 
associative ionization (AI) respectively. An example of (VIb. la) 
is 

He*(ls2s,^S) + Ar + "Koj He + Ar^(3p^,^P) + e (e) . (VIb. 2) 

There are basically two kinds of field-induced ioniza- 
tion processes : field-assisted (FA) and field-modified (FM) . 

FA processes occur in cases where the electronic continuum would 
not have been energetically accessible in the absence of the 
field, whereas FM processes occur in cases where field-free 
ionization is already possible, the field serving to bring new 
regions of the electronic continuum into accessibility. (VIb. 2) 
is an example of FM ionization. The system can be described by 
two potential curves: V^j for He*+Ar and V^ for He+Ar , V^j is 

embedded in the electronic continuum of V.j. and hence field-free 
ionization is possible. The introduction of the field enables 
new regions of the continuum to be accessible by collision, as 
illustrated in Figures 11 and 12. From these figures the qualita- 
tive picutre emerges that there will be three peaks in the 
emitted-electron energy spectrimi of the cross-section. The 
middle one is the field-free one, whereas the other two are 
fie Id- induced, spear ated from each other by 2iro). 

In the field-free case the bound-continuum aspect of 
the problem consists of the fact that the discrete state j<P(a>, 
corresponding to the surface V,^, can interact with the continuum 
of states j<j> .(.,£> corresponding to the surfaces V^+e, where £ is 
an arbitrary electronic energy. Any discretization procedure 
amounts to replacing the infinite set of continuum states |<}>^,e> 




He'(ls2s,^S)*Ar+-ftcu 


- He’(ls2s^S)*Ar 


21 - 


Fig. 11. Relevant potential energy curves for He* + Ar + "Ro) 
collisions . Emission of an electron with energy £+ requires 
absprption of a photon. 


by an infinite set of discrete states such that the 

Coulomb couplings, and V^g, leading to ionization are related 
by 

.00 

de a-.(e)V- , (VIb.3) 

3 £ de 

where a~(e) = <^ , eU , s> . (VIb.4) 

£ + + 



ag(£)' must satisfy the orthonormality condition 

de a*(e)a (e) = , 

Jq £ £ ££ 


(VIb.5) 





•He (ls2s, S)*Ar 




— He * Ar*(3p^,2p) 


Fig, 12. Relevant potential energy curves for_He* + Ar + "Rej 
collisions . Emission of an electron with energy £_ leads to 
emission of a photon. 


A particular discretization procedure is then specified by a par- 
ticular choice for the set a„{e). Very often the simple choice, 
satisfying (VIb.5) approximately, is used: 

a., (e) = (Ae) , for £ - 4 A&$&<£ + -^ As ; 

= 0 . , otherwise , (VIb.6) 

where the stepsize Ae can be chosen to be arbitrarily small. If 
one firrther introduces the Franck-Condon approximation which im- 
plies that only localized bound-discretized state couplings are 
significant, the diabatic representation of the potential sur- 
faces for nuclear motion with angular momentum L (11=1) is simpli- 
fied to 




(VIb.7) 



and y is the reduced mass of the collision system. The coupling 
Vjjg induces an avoided crossing where the curves and te 

cross, and the S-matrix describing transition to a discretized 
state, , can be computed by the semiclassical formalism 

described earlier. The S-matrix for transition to the true con- 
tinuum is then given by 

= y a„(£)Sp^ . (VIb.lO) 

£ ^ £ £ 

E 

When a field is present, two new regions of the con- 
tinuum are accessible. The first one is such that cross 

(N-D'Ko), and the second one such that cross 

e+(N+l)'R{». The new regions are such that both crossings occior 
where the (field-free) Coulomb coupling is dominated by the 
longer-ranged radiative coupling. Ionization in the first region 
requires absorption of a photon while ionization in the second 
region leads to emission of a photon. Semiclassical treatments 
can be applied to these field- induced cases entirely analogously 
to the field-free case. 

Vic. Time-Dependent Perturbation Theory 

Let us consider the following case of collisional exci- 
tation in the presence of a laser, 

A* + B + -Klo ^ A + B** , (VIc.l) 

where the asymptotic energy levels of the atoms A and B are il- 
lustrated in Figure 13. The formulation of the problem given 
below -is due to Weiner. 




Pig, 13. Asymptotic energy levels of a diatomic collision. The 
dotted line represents a virtual atomic state of B. 


A long-range dipole-dipole interaction between A and B 
produces a virtual B state at energy equal to the energy differ- 
ence between A and A*, which is a linear combination of the states 
corresponding to B and B* (mixing s and p characters) . Radiative 
coupling then carries the system from the virtual state to the 
final state B**. Direct radiative coupling (dipole coupling) 
between B and B** is symmetry forbidden (since both are S states) . 
Both the dipole-dipole and the radiative dipole interactions are 
assumed to be weak so that perturbation 'theory can be applied. 
These interactions are respectively given by 
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12 12 


12 


(Vic . 2)' 


23 


23 

Eh coscot 

O 


(Vic. 3) 


where E is the field strength of the laser and the u are the 
dipole moments of the respective atoms: 


= <A(s) [ vIa* ( p)> , 

A 


(Vie. 4a) 


= <B(s) |'u|b*(p)> , 

13 


(Vie. 4b) 


= <B* (p) (s) > ; 

5 


(vie. 4c) 


and R is the internuclear distance between A and B. The coupled 
equations (through second order) for the time-dependent expansion 
coefficients of the atomic wave function for B are (b^*^^=0/ 

bf)=0): . " 


= 0 

.(1) _ 1 , (0) /“l2^ 

^2 - iff ^ "^12^ 


(Vic. 5) 


(Vie. 5) 


,(2) ^ _L /“23^ 

3 iff ^2 ^23 


(Vic. 7) 


Putting b{°^=l and using the straight line approximation 
R(t)=vt with impact parameter b, we have 


( 1 ) 


12 12 
^A ^B 
iff 


^“12"^ 


dt 


/■u2_l 2^2, 3/2 
(b +v t ) 


12 12 
2 l^A^B 


iff 2, 

V b 


"U*"! 


“ 12 *" 


asymptotically 


12 12 
^^A ^B 


iiT 


vb 


(Vic. 8) 



where Kj_ is a first-order modified Bessel function and the 
asymptotic limit holds when Substituting (Vic. 8) into 

(Vic. 7) and using (Vic. 3), we have 
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•h 


„ 12 12 23 

% 

vb^ 


iAoit 


(Vic. 9) 


where the RWA has been applied and Hence 

„ 12 12 23 

(2) 1 ^A iAcot 

b ^ — (e -1) , 

±fi (Aw)vb 

and the finite-time transition probability is given by 
P(t) 


(Vic. 10) 




fl l^A 1 


ro 

•KAu J 


(1-cosAtut) . 

(Vic. 11) 
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P (t) is proportional to E , i.e./ linear with respect to the in- 
tensity of the radiation Chnplying a single-photon process) . 



PROBLEMS 


1. Using (Ib.4) to (Ib-8) verify (lb. 9 ) , in which the matrix 
elements are given by (Ib.lO) to (Ib.l2) . 

2 . Show that the diagonal elements of P [defined in (Ib.15)] 
vanish identically regardless of representation and that 

P is hermitian . 

3. Verify (Id.24) . 

4. Give an approximate, perturbative derivation of the Landau- 
Zener formula [(IIa.8)] in the low velocity limit. 

5. Following Stuckelberg ' s procedure write down the fourth- 

order differential equation for [where X-j_ and X 2 satisfy 
(IIc.l)] . Then assume a general WKB solution of the form 
(lie. 5) and obtain equations for S^ and S^ when coefficients 
of the powers of -fi are separately equated to zero. Finally, 
retaining only S and S^ in the expansion (lie. 5) derive the 
general WKB solutions for X-j_ X 2 • 

6 . In the Meyer-Miller treatment establish the equivalence 
between Hamilton's equations of motion for the classical 
electronic degrees of freedom [(Ile.lOb), (lie. 10c) and 
(Ile.lOd)] and the time-dependent coupled equations (Ile.l) 
with the trajectory R(t) determined by Ehrenfest's theorem. 

7 . Show that for the three-state nonadiabatic transition pro- 
blem the branch point structure of the adiabatic surfaces 
is still derived from square root functions, 

8. Show that C (R) [given by (Ib.25)] with a given by (llld.lO) 
diagonalizes U'^(R) as given by (IIId.7) . 

9. For a simple harmonic potential verify (IVa.7) and (IVa.8) . 

10. Verify (IVb.l2), (IVb.l3) and (IVb,14) . 

11. Write down an expression for the nonadiabatic transition 
S-matrix for the case of a one-dimensional two-state, 
two-intersection-point system. 

12. Derive (Va.9) . 

Using the Miller-George theory in the context of the two- 
state electronic-field representation, show that in the 
weak field limit the local transition probability is propor- 
tional to the field intensity. 


13 . 



14. Assuming ■ single-mode emission write down the matrix for the 
diahatic representation of the emission-electronic-field 
surfaces and simplify. 

15. Verify (VIb.3) and (VIb.5) .. 
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